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A  differential  geometric  description  of  crystals  with  continuous  distributions  of  lattice  defects  and  undergoing  potentially 
large  deformations  is  presented.  This  description  is  specialized  to  describe  discrete  defects,  i.e.,  singular  defect  distribu¬ 
tions.  Three  isolated  defects  are  considered  in  detail:  the  screw  dislocation,  the  wedge  disclination,  and  the  point  defect. 
New  analytical  solutions  are  obtained  for  elastic  fields  of  these  defects  in  isotropic  solids  of  finite  extent,  whereby  terms 
up  to  second  order  in  shain,  involving  elastic  constants  up  to  third  order,  are  retained  in  the  stress  components.  The  strain 
measure  used  in  the  nonlinear  elastic  potential  -  a  symmetric  function,  expressed  in  material  coordinates,  of  the  inverse 
deformation  gradient  -  differs  from  that  used  in  previous  solutions  for  crystal  defects,  and  is  thought  to  provide  a  more  re¬ 
alistic  depiction  of  mechanics  of  large  deformation  than  previous  theory  involving  third-order  Lagrangian  elastic  constants 
and  the  Green  strain  tensor.  For  the  screw  dislocation  and  wedge  disclination,  effects  of  core  pressure  and/or  possible  con¬ 
traction  along  the  defect  line  are  considered,  and  radial  displacement  contributions  arise  that  are  absent  in  the  linear  elastic 
solution,  affecting  dilatation.  Stress  components  are  shown  to  differ  from  those  of  linear  elastic  solutions  near  defect  cores. 
Volume  change  from  point  defects  is  strongly  affected  by  elastic  nonlinearity. 
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1  Introduction 

Analysis  of  defects  in  deformable  crystals  has  been  a  principal  application  of  differential  geometric  methods  to  continuum 
mechanics  of  solids  since  the  mid  twentieth  century.  The  earliest  works  linking  continuous  distributions  of  dislocations  to 
geometric  concepts  are  often  attributed  to  Kondo,  Bilby,  Kroner,  and  co-workers  [8,49,51].  Accompanying  such  works  were 
development  of  a  multiplicative  decomposition  of  the  deformation  gradient  [10,51]  into  lattice  deformation,  perhaps  more 
often  termed  “elastic”  deformation,  and  that  due  to  dislocation  glide,  perhaps  more  often  termed  “plastic”  deformation, 
and  a  linear  connection  with  vanishing  curvature  whose  torsion  can  be  related  to  the  dislocation  density  tensor  [8,49,51]. 
These  theories  established  a  mathematical  foundation  for  kinematics  of  finite  deformation  plasticity  theory  developed  in 
subsequent  decades  [6, 37, 55, 75, 80]. 

Early  geometric  continuum  theories  admitting  distributions  of  more  general  kinds  of  defects,  including  disclinations 
and  point  defects  such  as  vacancies,  substitutional  atoms,  or  interstitial  atoms  [4,5,9,51,60]  were  introduced  soon  after 
dislocation  theories.  In  some  such  descriptions  [9,50,51]  additional  terms  are  incorporated  multiplicatively  into  the  de¬ 
formation  gradient  to  represent  contributions  of  particular  mechanisms  distinct  from  dislocation  slip  (e.g.,  dilatation  from 
point  defects).  The  linear  connection  describing  geometry  of  the  defective  crystal  lattice  may  have  curvature  that  can  be  re¬ 
lated  to  disclination  density  and  non-metricity  that  can  be  related  to  inhomogeneous  distributions  of  point  defects.  General 
and  rigorous  geometric  analyses  of  heterogeneous  bodies  undergoing  finite  incompatible  deformation  appeared  around  the 
same  time  [12,64,88]. 

More  recent  works  have  sought  to  extend  or  even  reformulate  continuum  field  theories  of  defective  crystals  applicable 
to  the  finite  strain  regime.  Several  works  particularly  relevant  to  developments  in  this  paper  are  mentioned  in  what  follows; 
a  review  of  historic  and  more  recent  literature  is  beyond  the  present  scope  but  can  be  found  in  [77].  In  [I],  an  elastic- 
plastic  dislocation  field  theory  was  developed  that  reduces  to  nonlinear  elastic  theory  of  distributed  dislocations  [91]  in 
the  time-independent  limit.  Gauge  field  theory  has  been  applied  to  describe  lattice  defects  [34,53,58],  including  single 
dislocations  and  dislocation  distributions.  Dislocations,  disclinations,  and  point  defects  have  recently  been  analyzed  in 
theoretical  frameworks  incorporating  Cartan’s  moving  frames  [92-94].  Continuum  dislocation  theory  has  also  been  applied 
to  shells  undergoing  finite  deformation  [3,96]. 
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In  [19],  a  kinematic  description  was  developed  that  accounts  for  simultaneous  existence  of  dislocations,  disclinations, 
and  point  defects.  This  treatment  extended  previous  theory  [5,51,60]  by  including  statistically  stored  defects  [i.e.,  disloca¬ 
tions  (disclinations)  with  no  net  Burgers  (Frank)  vector]  and  incorporating  a  third  term  in  the  multiplicative  decomposition 
associated  with  residual  lattice  deformation  remaining  in  an  element  of  material  upon  stress  release,  induced  by  microscopic 
elastic  fields  of  defects  within.  Further  developments  on  the  origin  of  this  third  term  in  dislocated  or  twinned  single  crystals 
and  polycrystals  are  reported  in  [13, 17,22-24].  Thermodynamics  and  balance  laws  for  a  micropolar  medium  with  disloca¬ 
tions  and  disclinations  were  developed  in  [25],  and  the  geometric  treatment  was  extended  to  dielectric  solids  with  mobile 
vacancies  and  disclinations  associated  with  domain  walls  [21].  Further  derivations  of  fundamental  geometric  relations  and 
potential  applications  in  modern  crystal  plasticity  theory  are  reported  in  the  monograph  [14]. 

In  this  paper  a  general  theory  for  continuous  defect  distributions  is  presented  first  and  then  specialized  to  isolated  defects 
in  the  Eulerian  setting,  differing  from  works  by  Kroner  and  Seeger  [52]  and  Willis  [91]  that  dealt  only  with  dislocations  in 
the  Lagrangian  setting.  For  elastic  bodies  containing  dislocations,  disclinations,  and  point  defects,  some  parallels  between 
continuous  and  discrete  defect  theories  were  noted  by  De  Wit  [31,32],  albeit  restricted  to  linear  theory,  i.e.,  small  deforma¬ 
tion.  Nonlinear  theory  developed  in  the  present  work  extends  [31, 32]  by  associating  certain  singular  part(s)  of  total  finite 
distortion  with  inelastic  deformation  due  to  isolated  defects,  making  use  of  generalized  functions  [39]  to  address  singular 
quantities. 

Linear  elasticity  is  typically  deemed  sufficient  for  representing  elastic  fields  beyond  some  distance  from  the  core  of  the 
defect.  For  defect  densities  that  are  not  too  large,  the  mean  spacing  between  defects  is  great  enough  that  stress  fields  can 
be  addressed  using  linear  elastic  solutions.  However,  close  to  the  core,  strain  component(s)  from  linear  elastic  solutions  are 
quite  large,  rendering  the  linear  analysis  inaccurate  [38,81].  Furthermore,  some  physical  features  apparent  in  experiments 
and  atomic  simulation  such  as  residual  lattice  dilatation  [44,76,95]  and  phonon  scattering  [71]  cannot  be  predicted  by  linear 
theory,  but  can  be  predicted  using  nonlinear  elasticity. 

Analytical  solutions  to  nonlinear  boundary  value  problems  for  crystals  with  individual  defects  can  be  considered  rela¬ 
tively  scarce  relative  to  abundant  linear  elastic  solutions.  Exact  analytical  solutions  appear  to  exist  only  for  solids  described 
by  strain  energy  potentials  that  may  be  physically  unrealistic  for  real  crystals.  Exact  closed-form  solutions  can  be  obtained 
for  incompressible  elasticity,  as  has  been  demonstrated  for  the  screw  dislocation  [1,69,92,96].  It  has  been  noted  [38]  that 
displacement  fields  of  a  screw  dislocation  and  a  wedge  disclination  are  a  subset  of  a  larger  family  of  controllable  deforma¬ 
tions  for  which  universal  solutions  exist  in  incompressible  isotropic  elastostatics  [89].  Incompressible  constitutive  models, 
by  definition,  cannot  address  residual  elastic  dilatation  and  hence  are  not  physically  realistic  in  this  regard.  Analytical  so¬ 
lutions  for  wedge  disclinations  and  point  defects  have  been  reported  for  materials  obeying  incompressible  neo-Hookean 
elasticity  [93,94].  Eor  compressible  solids,  several  exact  analytical  solutions  for  dislocation  and  disclination  configurations 
have  been  derived  using  finite  strain  measures  [33,96,97],  but  the  semi-linear  energy  potential  used  in  these  analyses  may 
not  address  fully  nonlinear  behavior  in  typical  crystalline  metals  and  minerals.  In  particular,  this  semi-linear  potential  de¬ 
pends  on  only  two  elastic  constants  rather  than  the  usual  five  for  a  second-order  nonlinear  hyperelastic  and  isotropic  solid. 
An  exact  solution  has  also  been  derived  for  the  screw  dislocation  in  an  isotropic  solid  with  a  nonlinear  elastic  potential 
depending  on  the  small  strain  tensor  [46]. 

Approximate  analytical  solutions  have  been  found  for  solids  obeying  more  physically  realistic  strain  energy  potentials. 
In  such  approaches,  a  nonlinear  elastic  potential  is  invoked,  but  higher-order  products  of  displacement  gradients  entering 
the  stress  held  are  truncated  above  a  certain  order.  The  qualiher  “second-order”  is  used  here  to  denote  an  elasticity  solution 
wherein  stress  depends  quadratically  on  strain,  with  the  elastic  potential  cubic  in  strain  and  including  elastic  constants  of 
up  to  order  three.  A  “hrst-order”  solution  denotes  a  linear  stress-strain  response  and  a  quadratic  elastic  potential,  i.e.,  usual 
linear  elasticity. 

Kroner  and  colleagues  [52, 67]  instituted  a  geometric  stress  function  method  for  obtaining  second-order  accurate  solu¬ 
tions  for  isotropic  dislocated  bodies  obeying  Murnaghan’s  Eulerian  elastic  potential  [61].  Stress  helds  were  reported  for 
dislocations  in  inhnite  bodies  and  in  hnite  bodies  with  stress-free  cores,  for  singular  densities  of  edge  and  screw  dislo¬ 
cations;  this  solution  technique  explicitly  yields  stresses  but  not  displacements.  Seeger  and  Mann  [72]  derived  a  solution 
for  the  screw  dislocation  with  a  hollow  (stress-free)  core  using  an  isotropic  second-order  Lagrangian  analysis,  i.e.,  strain 
energy  potential  in  terms  of  Green  strain  with  up  to  third-order  elastic  constants.  Willis  [91]  formulated  an  iterative  pro¬ 
cedure  involving  Green’s  functions  for  inhnite  anisotropic  bodies  with  dislocations,  but  did  not  address  core  boundary 
conditions.  An  approximate  solution  for  an  edge  dislocation  in  an  inhnite  medium  obeying  isotropic  Lagrangian  elasticity 
was  reported  [48],  following  general  techniques  for  analysis  of  planar  second-order  nonlinear  elastic  problems  [2].  Teodo- 
siu  and  Soos  [81,82]  formalized  a  general  iterative  procedure*  for  analysis  of  dislocations  in  nonlinear  elasticity,  similar  to 
those  implemented  in  [38,72,73,91],  demonstrating  certain  uniqueness  properties  of  hrst-  and  second-order  solutions.  A 
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Such  iterative  methods  for  solution  of  nonlinear  elastic  problems  were  apparently  first  reported  by  Signorini  [74]  and  were  applied  intuitively  to  a 
number  of  problems  by  Murnaghan  [62]. 
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second-order  solution  for  the  isolated  screw  dislocation  in  a  finite  cylinder  with  core  pressure  was  derived  in  [81],  incorpo¬ 
rating  Eulerian  coordinates  and  the  Lagrangian  energy  potential  mentioned  above.  Advantages  of  Eulerian  coordinates  over 
Lagrangian  coordinates  have  been  discussed  [83].  More  recently,  gauge  theory  was  used  to  obtain  a  second-order  solution 
for  a  screw  dislocation  in  isotropic  Eulerian  elasticity  [58]. 

For  an  isolated  wedge  disclination,  a  solution  incorporating  both  Eulerian  coordinates  and  Murnaghan’s  Eulerian  poten¬ 
tial  was  derived  [68, 86],  but  core  pressure  was  not  addressed.  For  an  isolated  point  defect,  Seeger  and  Mann  [72]  extended 
the  approximate  second-order  Lagrangian  analysis  of  Murnaghan  [62],  who  treated  a  thick  spherical  shell  with  pressurized 
interior  and  exterior  boundaries.  In  this  work  [72],  both  the  defect  and  medium  were  treated  as  elastic  bodies  with  po¬ 
tentially  different  second-  and  third-order  elastic  constants;  see  also  [57]  that  addresses  related  nonlinear  elastic  problems 
including  interactions  among  spherical  inclusions. 

As  inferred  from  above,  second-order  solutions  involve  elastic  constants  of  orders  two  and  three.  These  approximate 
solutions  are  usually  justified  by  the  limitation  that  complete  elastic  constants  of  order  four  have  rarely  been  available 
for  real  materials,  and  solutions  of  higher  order  would  require  knowledge  of  fourth-order  constants.  Possible  relationships 
between  fourth-order  constants  and  geometric  features  of  the  edge  dislocation  core  have  been  noted  [47].  Longitudinal 
fourth-order  constants  can  be  inferred  from  shock  compression  data  [16],  and  Density  Functional  Theory  (DFT)  can  used  to 
predict  fourth-order  constants  [63],  though  accuracy  of  such  predictions  is  uncertain.  Regardless,  when  strains  become  large 
enough  that  third-order  elastic  constants  are  insufficient,  it  becomes  prudent  to  study  core  fields  using  atomic  simulation 
[20,  26,  27,  76]  rather  than  traditional  continuum  elasticity.  Nonlocal,  strain  gradient,  and/or  gauge  theories  [53,  58]  of 
continua  may  also  enable  a  physically  realistic  description  (e.g.,  smoothed  or  non-singular)  of  defect  cores.  Second-order 
solutions  derived  later  in  this  work  are  intended  for  regimes  and  materials  where  a  third-order  elastic  potential  is  accurate 
and  exclude  the  (singular)  defect  core.  Higher-order  gradients  of  deformation  (e.g.,  as  treated  in  a  general  constitutive 
formalism  in  [79])  are  not  incorporated  in  the  elastic  potential  used  here. 

Previous  second-order  elastic  analyses  of  lattice  defects  mentioned  above  have  either  used  a  Green  Lagrangian  elastic 
potential  W (£!)  [14,62,87]  or  Murnaghan’s  Eulerian  potential  W (e)  [61],  the  latter  restricted  to  isotropic  solids.  [Precisely, 
letting  F'  denote  the  two-point  deformation  gradient  tensor.  Green  strain  is  £J  =  ^{F'^F  —  1)  and  Almansi  strain  is 
e  =  |(1  —  F^^ F^^).]  In  contrast,  new  nonlinear  elastic  solutions  derived  in  the  present  work  incorporate  the  less 
often  encountered  Eulerian  strain  energy  potential  W{D),  where  D  —  |(1  —  F~^F^^).  Symmetric  tensor  D  is  termed 
an  “Eulerian”  strain  in  material  coordinates  [28],  is  invariant  under  rigid  rotations  of  the  spatial  coordinate  frame,  and 
was  introduced  by  Thomsen  [85]  and  soon  thereafter  Davies  [29]  for  describing  the  nonlinear  thermoelastic  response  of 
pressurized  cubic  crystals.  Unlike  W,  potential  W  can  be  applied  to  anisotropic  solids  [66,90].  This  model  has  demonstrated 
greater  intrinsic  stability  [43]  at  large  shear  and  compression  than  usual  Lagrangian  elasticity  incorporating  W  that  tends 
towards  instability  at  large  uniaxial  compression  [16].  It  is  also  more  accurate  for  addressing  the  response  of  most  crystalline 
solids  under  hydrostatic  compression,  whereby  at  second  order  it  degenerates  to  the  Birch-Murnaghan  equation-of-state 
[1 1, 84]  for  cubic  crystals  and  isotropic  media.  In  contrast,  second-order  Lagrangian  theory  tends  to  underestimate  pressure 
at  large  compression  [45].  Improved  accuracy  of  W  over  W  for  describing  hydrostatic  compression  of  lower  symmetry 
crystals  has  been  reported  [66,90].  The  nonlinear  elastic  response  of  diamond  under  finite  strain  predicted  by  DFT  [63] 
is  also  better  represented  by  Eulerian  elasticity  incorporating  W  than  Lagrangian  elasticity  incorporating  W,  with  elastic 
constants  of  up  to  order  four  considered  in  each  representation.  Recent  analysis  [16]  has  empirically  demonstrated  other 
advantages  of  Eulerian  theory  (based  on  W)  over  Lagrangian  theory  for  homogeneous  static  compression  and  homogeneous 
shear  of  representative  cubic  or  isotropic  solids  with  Cauchy  symmetry  of  third-order  elastic  constants.  For  crystals  of  lower 
symmetry  (e.g.,  quartz  and  sapphire)  subjected  to  uniaxial  shock  compression,  advantages  of  Eulerian  over  Lagrangian 
theory  were  less  evident.  Solutions  derived  later  in  the  present  paper  are  the  first  known  for  dislocations,  disclinations,  and 
point  defects  incorporating  elastic  potential  W (D). 

Eulerian  potential  W{D)  is  used  to  obtain  new  second-order  analytical  solutions  to  three  discrete  defect  problems  in 
the  present  paper,  which  is  structured  as  follows.  In  Sect.  2,  continuum  defect  theory  is  concisely  reviewed,  refining  certain 
aspects  of  earlier  work  [14, 19].  In  Sect.  3,  the  screw  dislocation  is  analyzed.  In  Sect.  4,  the  wedge  disclination  is  analyzed. 
In  Sect.  5,  the  point  defect  is  analyzed.  For  each  kind  of  defect,  it  is  shown  how  continuum  field  theory  with  a  multiplicative 
decomposition  of  the  deformation  gradient  and  additive  decomposition  of  a  linear  connection  can  be  applied  to  the  limiting 
case  of  a  singular  defect  density.  A  second-order  elastic  analysis  is  used  in  each  case,  with  solutions  expressed  in  Eulerian 
coordinates.  For  the  dislocation  and  disclination,  core  pressure  is  included.  All  three  solutions  are  believed  to  be  new. 
Representative  calculations  demonstrate  differences  between  nonlinear  and  linear  elastic  solutions.  Conclusions  follow  in 
Sect.  6.  Appendices  A  and  B  provide  supporting  discussion  on  finite  deformation  in  cylindrical  and  spherical  coordinate 
systems  used  in  the  derivations,  while  Appendix  C  describes  nonlinear  elastic  potentials  and  associated  second-  and  third- 
order  elastic  constants  and  clarifies  differences  among  several  nonlinear  elastic  theories.  Direct  tensor  notation  [64,78]  and 
coordinate  notation  of  tensor  calculus  [36, 70, 77]  are  used  interchangeably  as  needed. 
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2  Continuous  defect  distributions 

2.1  General  continuum  theory 

Consider  a  solid  body  whose  unstressed  reference  configuration  Sq  is  perfectly  crystalline  and  simply  connected.  At  some 
time  t,  this  body  in  its  current  configuration  representation  05  may  contain  defects  and  may  support  stresses,  both  residual 
and  due  to  externally  applied  traction.  Deformed  body  55  is  assumed  to  remain  simply  connected  in  the  continuum  theory, 
and  defects  are  characterized  by  spatially  continuous  density  functions.  A  local  volume  element  is  associated  with  a  material 
particle  at  position  X  in  three-dimensional  Euclidean  space  and  is  mapped  to  position  x  at  time  t.  This  element  may  contain 
a  number  of  defects  such  as  dislocations,  the  product  of  whose  number  and  Burgers  vector  remains  finite  [91]. 

Henceforth  attention  is  restricted  to  a  deformed  configuration  at  a  single  time  t,  such  that  the  description  becomes  time- 
independent.  Let  {X^}  and  {x“}  denote  coordinate  charts  covering  55o  and  55,  respectively.  The  following  one-to-one  and 


presumably  three  times  differentiable  functions  exist: 

x  =  x{X),  X  =  X{x).  (2.1) 

Letting  9a(’)  =  and  da{-)  =  d{-)/dxa,  natural  basis  vectors  are 

GA{X)  =  dAX,  g,{x)  =  daX.  (2.2) 

The  deformation  gradient  and  its  inverse  are 

F  =  dAX(g>G^,  =  daX  (g)  g'^ ,  (2.3) 

where,  with  Gab  (-^)  and  gab(x)  referential  and  spatial  metric  tensor  components, 

Ga=GabG^  =  {Ga-Gb)G^,  ga  =  9ab9^  =  {ga  ■  gb)g'^ ■  (2.4) 

In  component  form, 

F^A=dAx\  F-\i=daX^.  (2.5) 

Conventional  gradient  operators  on  55o  and  55  are  denoted  by  Vo  and  V.  Lor  example,  letting  V  and  v  denote  material  and 
spatial  vector  fields,  covariant  derivatives  are 

VoV  =  dBVgG^  =  {ObV^  +rg^cy^')GA(gG^,  (2.6) 

Vv  =  dbv  gg^  =  (dbv‘‘  +  Jbcv‘')ga  0  g\  (2.7) 

where  Levi-Civita  connection  coefficients 

Tjgc  =  \G^^  {dBGcD  +  OcGbd  —  doGBc),  %'c  =  \g‘"^{,dhgcd  +  dcgbd  —  dagbc)-  (2.8) 


Lrom  the  properties  of  mappings  in  (2.1)  and  presumed  positivity  of  local  volume  elements  in  55o  and  55  in  right-handed 
coordinate  systems,  restrictions  det  >  0  and  det  >  Demerge.  Since  and  {a;“}  are  holonomic,  9^[£)b(-)]  = 

Ob  [94(’)]  and  c)o[i96(’)]  =  56['9a(’)]’  and  since  the  body  is  embedded  in  Euclidean  space,  Riemann-Christoffel  curvature 
tensors  constructed  from  (2.8)  vanish  identically.  Lrom  (2.5)  and  symmetry  of  connection  coefficients  in  (2.8),  the  following 
local  compatibility  (i.e.,  null  curl)  conditions  apply: 

VoxF  =  0,  Vx(F-i)  =  0.  (2.9) 

Lor  defective  crystals  a  multiplicative  decomposition  of  F  is  invoked  [10, 14,51]: 

F  =  F^F^,  F^a  =  “a  “a  •  (2- 10) 

Lattice  deformation  F^  accounts  for  elastic  stretch  (recoverable  and/or  residual)  and  rotation;  plastic  deformation  F^ 
accounts  for  slip  from  dislocations  that  preserves  the  lattice  structure  [9, 10].  Both  fields  are  assumed  invertible,  with  positive 
determinants,  and  at  least  twice  differentiable  with  respect  to  X^  and  a;“ .  Neither,  however,  is  necessarily  integrable  to  a 
motion  nor  satisfies  conditions  analogous  to  (2.9);  i.e.,  lattice  and  plastic  deformation  tensors  are  generally  anholonomic. 
Thus,  the  intermediate  configuration  18  of  the  body  implied  by  (2.10)  cannot  be  embedded  in  three-dimensional  Euclidean 
space  [15,41].  Letting  {ga}  denote  a  generally  non-coordinate  basis,  with  reciprocal  basis  {g“}, 

FL  ^  pL  F^  =  F^  ^_Aga  ®  G^ .  (2.11) 
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Introduced  at  each  point  X  G  $0  is  a  triad  of  undeformed  lattice  director  vectors  {d^}  and  their  reciprocals 
such  that  (d'"',  ds)  =  5g  with  (•,  •)  a  scalar  product.  Referential  director  vectors  are  orthonormal:  d,4  •  ds  =  Sab  and 
d.4  =  (^Asd^.  In  crystal  structures  with  cubic,  tetragonal,  or  orthorhombic  symmetry,  a  referential  lattice  vector  can  be 
assigned  parallel  to  each  edge  of  the  unit  cell.  Also  introduced  at  each  a;  €  S  is  a  triad  of  deformed  lattice  directors  {da} 
and  their  reciprocals  {d“},  such  that  (d“  ,  dj)  =  6^.  Deformed  director  vectors  are  generally  neither  orthogonal  nor  of  unit 
length: 


do  rriL  a  cq  jA  j  771L — la  cA  j 

=  F  ,  da=F  .aOadA- 

The  metric  tensor  associated  with  the  deformed  lattice  is 
d={da-di)g<^(S)g^  =  F^-^F^-\ 


In  components,  this  metric  and  its  inverse  are 

7  rpL — la  r  TTtL — 1/5  job  771L  a  socB  rriL  h 

dab  =  F  ,  d  =F  ^ F 


(2.12) 


(2.13) 


(2.14) 


leading  to 


da=da&d^  d“=d“'d6. 


(2.15) 


Components  of  lattice  deformation  can  be  inferred  from  scalar  products  of  referential  and  deformed  directors: 

F^  \  =  (d“  ,  dA)S^  ,  =  (da  ,dAS1, 

and  if  a  coordinate  frame  on  IB  is  chosen  parallel  to  the  reference  lattice  such  that  {ga}  — *■  {dAS^}, 

^  Qa®  d“ ,  =  da®  g‘". 


(2.16) 


(2.17) 


In  the  continuum  field  theory  of  defects,  the  geometry  of  the  (possibly  defective)  crystal  lattice  is  characterized  by 
covariant  derivatives  of  the  spatial  lattice  director  fields,  defined  as  [14, 19,60] 

Vida  =  d,,da  -  fl'^da,  V id!^  =  dfed"  +  f[,a“dC  (2.18) 

A  complete  description  of  the  deformed  lattice  is  given  by  the  fields  {d„}  and  {  Vda},  where  a  =  1,2,3.  From  (2.16) 
and  (2.18),  this  description  is  furnished  by  knowledge  of  the  fields  F^  (up  to  9  components)  and  F  (up  to  27  components) 
which  are  available  at  each  point  x  G  18,  or  equivalently  at  each  point  X  G  IBo  since  x  and  X  are  one-to-one  in  (2.1). 
An  additive  decomposition  of  connection  coefficients  associated  with  V,  assumed  at  least  once  differentiable  over  18,  is 
invoked  [14, 19,60]: 


-p..a  -pi. a  I  ^••o 

^bc  —  ^  be  V6c  • 


(2.19) 


Coefficients  f  are  spatial  representations  of  those  of  the  usual  non-Riemmannian  connection  of  dislocation  field  theory, 
referred  to  as  a  “crystal  connection”  [54]  in  the  context  of  nonlinear  elastoplasticity  or  an  example  of  a  “Weitzenbdek” 
connection  in  differential  geometry  [53,92]: 


=  -F^ 


(2.20) 


The  torsion  tensor  of  (2.19)  is 

T  =  f  ®g‘-'®ga,  fi-:  =  1  (f-  -  f  )  =  F^  J  -t-  , 


(2.21) 


where  indices  in  square  brackets  are  skew.  Verification  is  straightforward  that  the  Riemann-Christoffel  curvature  tensor 
from  (2.20)  vanishes  identically,  meaning  this  connection  is  flat  or  integrable: 


R  =  -f  ®g‘''®g‘^®ga 

=  2F^ ®g'''®g'^®ga  =  0.  (2.22) 
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Indices  in  vertical  bars  are  excluded  from  (anti)symmetry  operations.  The  curvature  tensor  of  (2.19)  is  [14,70] 

=  2(V[i(5"|^  +  T{,'cQ'ed  +  Q[b]d\Qcie)9'‘  ®  g^  ®  g'^  ®  ga, 


(2.23) 


where 


VbQTd  =  dbQTd  +  n!Qfd  - n!QTd  - ^TdQTe- 

Nonmetricity  of  connection  (2.19)  is 

^ cdab  ~  ^cdah  ^ ca  ^be  dae  —  {.Qca  ^eb  Qcb  *^ea)  ~  ‘^Qc{ab)  i 

where  indices  in  parentheses  are  symmetric  and  d  is  used  to  lower  indices  of  Q.  Covariant  components  of  R  are 

Rbcda  d^bcd  dea  ‘^{'^[bQcjda  ^bc  Qeda  Q[b\d\^c]aeb 

^[bc]{da)  ^hc{ad)  [bQ c]{da)  “f  Qe{da)b 

R[hc][da\  ^hc[ad\  [bQc][da]  “t”  ^bc  Qe[da])  Q6[(;Q|c|a]e  “t“  Qc[dQ\d\a]e- 


(2.24) 


(2.25) 


(2.26) 

(2.27) 

(2.28) 


Curvature  and  non-metricity  arise  from  Q  and  vanish  when  f  =  f ;  nonmetricity  further  depends  only  on  the  symmetric 
part  of  Q  in  (2.25). 

Mathematical  quantities  associated  with  defect  densities  are  now  defined.  The  spatial  dislocation  density  tensor  is 


a  =  -T-.e  =  t'^‘^^f-^fga®gb, 


where  e  is  the  spatial  permutation  tensor.  Usual  dimensionless  permutation  symbols  are  eabc  =  e“^'^  =  5“^/^ 
gil2^abc  g  _  (jet[c)aa;  •  dbx\.  Let  (•)^  denote  the  transpose  of  a  second-order  tensor.  The  spatial  disclination  density 
tensor  is 


(2.29) 


^abc  — 


0  =  -i(e:,R:e)T=  (2.30) 

Henceforth  restricting  nonmetricity  as 

Qa(bc)  =  -\^adbc  ='^adbc,  Qb(cd)d'^°'  =  \d°''^  {Qbcd  +  Qbdc)  = bd'^ ,  (2.31) 

a  vector  quantity  associated  with  point  defects  can  be  defined  [14, 21]: 

^  =  -VxT  =  e:  =  (2.32) 

Inverting  (2.29),  (2.30),  and  (2.32)  enables  reconstruction  of  the  torsion  and  curvature: 

Ti'c  =  \^<^bda°‘'^ ,  Rab[cd]=  (^bae(^cdfO^^,  Rab{cd]  =  ^bae'd'' dcd  ■  (2.33) 


Identify  18  as  a  manifold  equipped  with  connection  V.  Let  {u,v,w}  be  smooth  vector  fields  on  the  tangent  bundle  of 
<8.  Then 


(VVtt)  :  (v  (g)  lu  —  lu  (g)  u)  =  {Ru  —  2T'Vu)  :  {v  (g  lu),  (2.34) 

or  in  a  coordinate  chart  {a:“}  [70], 

(2.35) 

This  identity  can  be  used  to  prove  that  under  a  “Cartan  displacement”,  a  vector  u  undergoes  a  first-order  (with  respect 
to  transported  distance  da:)  change  upon  parallel  transport  about  an  infinitesimal  circuit  dc  of  the  form  m  — s-  it  +  5u, 
where  [60,70,89] 

5u  =  -{\  Ri"-  u<^  +  )da'"=  Qa .  (2.36) 
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The  oriented  area  element  enclosed  by  dc  has  magnitude  da,  unit  normal  n,  and  obeys 

da  =  da'"' 96  ®  Qc  =  e'’“nada  gt,  ®  Qc,  Uada  =  |ea6cda'"'.  (2.37) 

From  (2.36),  Su  consists  of  two  parts:  one  proportional  to  u  due  to  curvature  and  the  other  independent  of  u  due  to  torsion, 
the  latter  associated  with  “closure  failure”  of  the  circuit  itself.  Motivated  by  this  definition,  an  infinitesimal  spatial  Burgers 
vector  db  is,  with  a:  ^  a:  +  db  under  a  Cartan  displacement  [14, 60], 

db(a;)  =  -(f  +  ^Rx)  :  da.  (2.38) 

Extending  dc  to  a  finite  circuit  c  enclosing  area  a  (i.e.,  a  Burgers  circuit  on  the  deformed  crystal)  and  choosing  {ga}  constant 
over  IB,  the  total  spatial  Burgers  vector  is 


b=  /  db=  —  /(r+  kRx)  :  da 


This  can  be  expressed  in  terms  of  defect  densities  using  (2.33)  [14]: 


(2.39) 


b“  =  /  a“'’n6da+  /  nhda  +  I  x‘‘'d'’nbda 


(2.40) 


2.2  Particular  cases:  single  defect  types 

It  is  instructive  to  consider  cases  when  only  one  type  of  defect  is  present.  A  sufficient  description  of  a  crystal  containing 
only  dislocations  is  obtained  by  forcing  Q  =  0  in  (2.19),  such  thatF  — >  F  and  V  ^  V  describe  a  flat  connection  with 
non- vanishing  torsion  and  vanishing  nonmetricity.  Differentiation  of  (2.12)  with  cIa  =  constant  yields  the  parallel  transport 
equation 


di,da  =  =  f 6;^d,  O  Vftda  =  0.  (2.41) 

The  “local”  or  spatial  dislocation  density  tensor  becomes 

a  =  -F‘^[V  X  a"''  = 

Defining  the  “true”  dislocation  density  and  infinitesimal  Burgers  vector  as 
a  =  F^~^a  =  —V  X  db  =  dnda, 

application  of  Stokes’s  theorem  with  (2.1 1)  produces  the  true  Burgers  vector 

b  =  b"9a  =  f  db  =  j  anda  =  —  ®  F^~^dx  =  —  ®  F^dX,  (2.44) 

J  a  J  a  J  c  J  (t 

where  £  is  the  image  of  c  on  IBq  and  {ga  }  are  chosen  constant  over  the  domain.  For  a  density  of  dislocations  all  having  the 
same  unit  tangent  line  ^  and  discrete  Burgers  vector  b,  (2.42)  yields  a  =  pb  0  where  p  is  the  number  of  dislocations  per 
unit  area  a.  Accurate  when  p  is  not  too  large,  in  dislocation  theory  it  is  conventionally  assumed 

=  F^F^  «  F  =  F^F^F^  «  F^F^,  (2.45) 

where  F^  accounts  for  rigid  rotation  and  stretch  conjugate  to  applied  stress  [6,  55,  80],  and  F^  accounts  for  residual 
distortion  (e.g.,  dilatation  from  dislocation  cores)  [22]  not  addressed  by  slip  deformation  F^  which  is  by  definition  lattice 
preserving  [9, 10, 14,22]  and  thus  isochoric. 

For  a  body  containing  only  disclinations  with  no  history  of  slip  or  point  defects,  the  following  conditions  are  imposed: 

F  =  F^  =  F^F\  a  =  0,  =  0,  (2.46) 

where  F^  accounts  for  residual  deformation  associated  with  disclinations.  These  lead  to  the  following  constraints  that 
reduce  Q  to  no  more  than  6  independent  components: 

T;r=QiH=0,  Qaibc)=0-  (2.47) 


(2.42) 

(2.43) 
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Thus  f  and  V  now  describe  a  Riemannian  connection,  i.e.,  vanishing  torsion  and  nonmetricity  but  non-vanishing  curvature. 
Thus  a  metric  exists  whose  Levi-Civita  connection  is  F;  this  metric  is  an  integral  function  of  Q  that  degenerates  to  the  Levi- 
Civita  connection  of  dab  when  Q  vanishes.  The  local  Burgers  vector  is 

S=  /  :  (0n(S)  x)da.  (2.48) 

J  a 

For  a  density  of  disclinations  all  having  the  same  unit  tangent  line  ^  and  discrete  Frank  rotation  vector  u,  the  spatial 
disclination  density  tensor  is  0  ~  pu)  ^  where  p  is  the  number  of  disclinations  per  unit  area  a.  The  inhnitesimal  Burgers 
vector  then  reduces  to  db  =  pd~^  (a)  x  x)  (^,  n)da.  Unlike  the  case  for  dislocations  only,  here  the  Burgers  vector  depends 
on  choice  of  coordinate  system  by  which  x  is  measured,  i.e.,  db  is  not  translationally  invariant.  According  to  (2.48),  each 
disclination  is  effectively  assigned  a  radius  [14,25,60]  measured  from  the  origin  of  this  coordinate  system. 

For  a  body  containing  only  point  defects,  with  no  history  of  slip  or  disclinations,  the  following  conditions  are  imposed: 

F  =  =  F^F\  a  =  0,  0  =  0,  (2.49) 

where  F^  accounts  for  residual  deformation  associated  with  point  defects,  which  is  hereafter  assumed  dilatational  or  spher¬ 
ical: 


=  +  =  (2.50) 

The  number  of  point  defects  per  unit  reference  volume  is  p,  and  6v  is  the  residual  volume  change  per  defect.  The  number 
of  defects  per  unit  intermediate  volume  is  p,  and  p  =  p/(l  +  p5v).  The  following  constraints,  necessary^  for  satisfaction 
of  (2.49),  are  imposed,  corresponding  to  integrable  F^~^  and  vanishing  (2.21)  and  (2.28): 

r,”"  =  0  ^  -  Qcdt  -  Qbcd  +  Qbdc). 

A[bc][da]  =  0  2V [bQc][da]  —  ‘^Q\^bc]Qe[da]  H“  Qb[dQ\c\a]e  Qc[dQ\b\a]e’  (2-51) 

The  local  Burgers  vector  is 

b=  /  x{'0,n)da.  (2.52) 

J  a 

As  was  the  case  for  disclinations,  here  the  Burgers  vector  is  not  translationally  invariant.  The  Burgers  vector  from  (2.52)  is 
proportional  to  distance  from  the  origin  and  is  directed  radially  outward  from  the  origin.  If  the  density  of  point  defects  is 
uniform  throughout  an  inhnite  body  03,  then  (2.52)  must  vanish  for  any  choice  of  coordinate  system,  which  implies  r?  =  0 
for  any  such  uniform  distribution. 

In  general,  when  dislocations,  disclinations,  and  point  defects  are  all  present,  three-term  decomposition  F  =  F^F^F^  = 
F^F^  applies.  The  primary  distinction  between  F^  and  F^  is  that  the  former  affects  the  lattice  directors  through  (2.12) 
while  the  latter  does  not.  Inelastic  deformation  associated  with  distributed  disclinations  and  point  defects  includes  di¬ 
latation  (as  well  as  possible  shape  change  and  rotation)  and  cannot  be  lattice-preserving;  therefore,  such  deformation  is 
encompassed  by  F^.  Inelastic  deformation  due  to  dislocation  glide,  i.e.,  shearing  on  slip  planes,  is  both  isochoric  and  lattice¬ 
preserving,  and  is  measured  by  F^.  However,  in  geometrically  linear  theories  of  crystal  plasticity,  wherein  slip  directors 
are  assumed  constant,  inelastic  deformation  from  disclinations  and  point  defects  is  often  assigned  to  plastic  distortion  along 
with  deformation  from  slip  [14,31,32]. 

2.3  Nonlinear  elasticity 

Applications  of  the  theory  in  subsequent  sections  consider  isolated  line  defects  (single  dislocation  or  disclination)  and 
isolated  point  defects  (single  vacancy  or  substitutional/interstitial  atom).  The  method  of  solution  involves  assignment  of  a 
singular  density  of  defects  at  a  line  or  a  point  and  then  excluding  a  small  region  (e.g.,  a  cylinder  surrounding  the  line  or 
a  sphere  surrounding  the  point)  from  the  domain  over  which  equilibrium  equations  are  solved.  For  such  applications^,  the 

^  Previously  [14, 19]  conditions  Qa[(,c]  =  0  were  imposed  for  a  body  both  point  defects  and  (possibly)  dislocations;  these  conditions  are  sufficient 
but  not  necessary  for  vanishing  disclination  density. 

^  In  prior  work  [14,  25]  considering  time  dependent  problems  involving  continuous  defect  distributions,  dependence  of  strain  energy  on  defect 
densities  was  permitted,  and  additional  balance  laws  for  thermodynamic  conjugates  (e.g.,  couple  stresses)  to  defect  densities  were  proposed.  These 
additional  balance  laws  are  trivially  satisfied  in  the  present  context  since  defect  densities  vanish  over  domains  excluding  defect  cores,  but  in  other 
problems  they  are  intended  to  affect  kinetics  and  restrict  equilibrium  defect  distributions.  Proper  accounting  for  stored  energy  of  cold  work  [7]  in 
an  externally  unstressed  element  of  material  requires  that  the  element’s  strain  energy  depend  on  defect  density  [14]. 
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following  hyperelastic  constitutive  model  for  strain  energy  per  unit  reference  volume  W  suffices: 

W  =  W(F^),  P  =  dW/dF^  =  Ffg'^  0ga.  (2.53) 

Here  P  is  a  kind  of  (generally  non-symmetric)  first  Piola-Kirchhoff  stress.  To  ensure  invariance  under  rigid  rotation,  it  is 
implicitly  understood  that  W  depends  only  on  some  measure  of  stretch  associated  with  F^.  The  usual  symmetric  Cauchy 
stress  and  corresponding  static  equilibrium  equations  on  S,  in  the  absence  of  body  force,  are 

a  =  J-^PF^'^,  Va^O,  (2.54) 

where  the  Jacobian  determinant  J  =  ~  Equilibrium  equations  can  be  mapped 

to  the  intermediate  configuration  18  as  demonstrated  in  coordinate-free  notation  in  [59, 64]  and  for  particular  choices  of 
anholonomic  coordinates  in  [15];  torsion  of  the  crystal  connection  and  possible  skew-symmetric  parts  of  the  coordinate 
connection  on  B  enter  the  result. 

3  Screw  dislocation 

3.1  Geometric  description 

Consider  a  deformed  body  18  with  boundary  018.  In  the  reference  configuration,  the  image  of  this  body  and  its  boundary 
are  18o  and  018o.  Cylindrical  coordinates  covering  18  and  18o  are  {r,9,z)  and  {R,Q,Z)  and  are  described  fully  in  Ap¬ 
pendix  A  (See  Supplementary  Material,  Wiley  online  library),  including  basis  vectors;  Cartesian  coordinates  are  {x,  y,  z) 
and  {X,Y,Z). 

An  Eulerian  description  of  deformation  is  invoked  here,  essentially  inverting  the  Lagrangian  description  of  Seeger  and 
Mann  [72].  Let  the  unit  tangent  line  of  the  dislocation  be  oriented  along  the  z-axis,  i.e.,  ^  =  gz-  The  general  ansatz  for 
deformation  due  to  a  screw  dislocation  in  an  isotropic  body  is 

R(r)=r  —  u{r),  0(0)  =  0,  Z{6,z)  =  z  —  w{ip)  —  s{z).  (3.1) 

Radial  displacement  is  u{r).  Angle  6  is  multi-valued  upon  complete  revolutions  around  the  dislocation  line.  Slip  plane  © 
corresponds  to  the  half-plane  (y  =  0,  x  <  0),  or  0  =  ±7r.  Following  [32],  ip  is  introduced  to  maintain  uniqueness  of  the 
deformation  field  and  enable  description  of  the  slip  discontinuity  via  generalized  functional  analysis  [39].  Angle  ip  is  single 


valued  but  discontinuous,  and  is  defined  as  [32] 

ip{x,y)  =  tan“^(?//x)  -f  7rH(-x)[H(y)  -  H(-y)],  (3.2) 

where  the  arctangent  is  restricted  to  its  principal  value  in  the  range  (—  f ,  f )  and 

H(x)  =0  Vx  <  0;  H(x)  =  1  Vx  >  0.  (3.3) 

Note  that 

y:  e  (-tt.tt),  |(/?]  =  V5|9=7r  -  </5|9=-7r  =  27r,  If  =  9  V0e(-7r,7r).  (3.4) 

Generalized  derivatives  [39]  of  (p(x,  y)  are  [32]: 

dx<^  =  -y/r‘^,  dyip  =  x/r"^ +  2T:n{-x)  6  {y),  (3.5) 

where  S{y)  —  dH(?/)/dy  =  is  the  Dirac  delta  function.  The  following  generalized  derivatives  are  also  noted  for 

functions  whose  ordinary  derivatives  are  0(l/r^)  [32]: 

dx{x/r^)  =  l/r"^  —  2x‘^ /r'^  +  Tr6{r),  dy{y/r‘^)  =  l/r"^  —  2y'^ /r'^  +  Tr5{r),  (3.6) 

where  S{r)  =  6{x)5{y)  in  polar  cylindrical  coordinates. 

For  the  screw  dislocation,  the  true  Burgers  vector  and  function  w  are 

b  =  bgx  =  bGz,  w  =  b(p/{2'K).  (3.7) 

The  true  Burgers  vector  is  constant.  Since  0  =  0,  basis  vectors  in  Appendix  A  (See  Supplementary  Material,  Wiley  online 
library)  here  obey  {l/r)gg  =  {1/R)Gq  =  eg  =  Eq.  Displacement  u  =  x  —  X  in  (A. 38)  reduces  to 

u  =  u{r)gr  +  [w(v?)  +  s{z)]gz  =  u{r)GR  +  +  s{z)]Gz-  (3.8) 
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Displacement  is  continuous  except  for  that  due  to  the  jump  in  w  across  6: 

|m]  =  =  b.  (3.9) 

In  curvilinear  coordinates,  the  inverse  deformation  gradient  from  (3.1)  is 

F-i  =  F-^^Gr  0  9’'  +  F-i®Ge  0  g®  +  Gz  0  Gz  0  g® 

=  (di?/dr)G/j  <8>  g""  +  G©  <8>  g^  +  d^ZGz  <8>  g^  +  dgZGz  <8>  g^ 

=  [1  -  u'(r)]Gij  (g)  g’'  +  Ge  (g)  g^  +  [1  -  s'(z)]Gz  gg""  -  dgwGz  gg^,  (3.10) 

where  it'  =  dii/dr  and  s'  =  ds/dz.  From  (3.5), 

F~^g  =  —  r  sin^cij;!!!  +  r  cos6*c)„it;]  =  —  [1  —  27rrH(— a;)i5(7/)].  (3.11) 

■  ZTT  ZTT 

To  fix  translational  invariance  of  the  origin  and  ensure  volume  remains  positive, 

m(0)  =  0,  s(0)  =  0;  (l-it/^)(l-it')(l-s')  >  0.  (3.12) 

The  final  decomposition  (2.45)  is  now  applied,  such  that 

F  =  F^F^,  F~^=F^-^F^-\  (3.13) 

Extending  the  linear  theory  of  De  Wit  [32]  to  the  geometrically  nonlinear  regime,  the  “elastic”  part  of  F~^  is  defined  as 
that  which  is  continuous  and  single- valued  over  IB  except  at  r  =  0: 

FE-1  =  [1  -  u'{r)]GR  g  g'-  +  Ge  g  g®  +  [1  -  s'{z)]Gz  g  g"  -  [V(2^)]Gz  g  g^  (3.14) 

It  follows  from  (3.1 1)  and  (3.13)  that 

=6rH(-tc)5(i/)Gzgg®  +  F®-i  =  [6rH(-:r)5(i/)Gz  g  (F®  V)  +  •  (3.15) 

In  the  reference  configuration,  the  image  of  the  slipped  surface  ©o  is  (W  <  0,  F  =  0),  on  which  r{R)  and  0(0)  are 
continuous,  but  Z{9,  z)  is  discontinuous.  Let  m  denote  a  normal  to  ©  directed  from  0  =  — tt  to  0  =  +7r: 

m{r,R)  =  —rg^-F^  =  —rG^=—{l/R)gg  =  —{r/R)eg  on  ©.  (3.16) 

The  plastic  part  of  (3.15)  is  the  simple  shear 

F^-\x)  =  l-H{-x)S{y)b(Zm,  F^{X)  =  1  +  R{-X)S{Y)b  (Z  m.  (3.17) 

Since  b  m  =  —{br/R)ez  -  eg  =0,  slip  is  isochoric:  det  F^  =  1.  Integration  of  F^  over  a  rectangular  volume  V  traversed 
by  ©0  produces,  with  the  approximation  r/R  «  1, 

(1/F)^  F^dV  ^l  +  V~^b(Zm  R{-X)5{Y)dXdYdZ  =  1  + -fs  ®  m,  (3.18) 

where  7  =  bA/V  is  shear  in  direction  s  =  g^  due  to  slip  over  area  A  =  f  li(—X)dXdZ.  From  the  above  description, 
F{X)  =  F^  VXg<Bo\©o,  F{X)  =  F^F^  on  ©0; 

F^\x)  =  F^-'^  Va:  G  *B\©,  F~\x)  =  F^-^ F^-'^  on  ©.  (3.19) 

It  becomes  prudent  to  select  coincident  cylindrical  coordinate  systems  on  IB  and  IBq,  such  that  g^  (X)  =  S^Ga  {X)  [15], 
and  in  these  coordinates 

FE  =  F®“„g,gg“  =FE“^g,  (gG'\  VX  G  95o\©o.  (3.20) 

The  dislocation  density  tensor  can  now  be  computed.  In  Cartesian  coordinates,  using  (3.5), 

F~\x  =  %/(27rr^),  F~\y  =  -b[x/{2Trr^)  +  H(-a;)(5(g)].  (3.21) 
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Decomposing  this  into  continuous  and  singular  parts, 

=  W(27rr2),  F^-^%  =  -bxl{27Tr^),  (3.22) 

The  only  non- vanishing  component  of  true  dislocation  density  is,  from  (2.43)  and  (3.6), 

^Zz  ^  ^  g^pE-iz^  _  Q^pE-iz^  ^  [b/(2n)][dy{y/r^)  +  d:,{x/r‘^)]  =  bd{r).  (3.23) 

Therefore, 

a  =  b6{r)Gz  <^gz  =  5{r)b(^^.  (3.24) 

Consider  a  closed  loop  c  encircling  the  dislocation  line;  area  a  enclosed  by  this  loop  has  constant  normal  vector  n  =  =  ^. 

Applying  Stokes’s  theorem,  consistency  with  (2.44)  with  is  demonstrated: 

b=  6{r)b{^,n)rdrd6  =  —  /  (V  x  F^“^)nda 


=  -  6  F^-^dx  =  -  (b  F 


:^E-1Z 


gGzd0  =  —  (p  de  =  b.  (3.25) 

27r  J, 

The  second  of  (2.9)  can  be  verified  by  taking  the  curl,  in  the  sense  of  a  generalized  derivative  [39],  of  (3.15): 

-V  X  (F~i)  =  -V  X  (F^-^)  -f  V  X  [H(-a:)(5(t/)b  O  m]  =  a  -  S{-x)S{y)b  O  =  0.  (3.26) 


3.2  Nonlinear  elastic  analysis  and  general  solution 

A  nonlinear  elastic  boundary  value  problem  for  the  screw  dislocation  is  constructed  as  follows.  Let  body  05  consist  of  a 
cylinder  of  outer  radius  Rq  with  the  dislocation  line  located  along  r  =  0.  A  cylindrical  core  region  r  <  ro  is  removed  from 
the  body  and  excluded  from  the  solution  domain.  The  length  of  the  cylinder  is  L,  such  that  the  body  is  contained  within  the 
limits  —  j  <^<  h  though  the  forthcoming  analysis  also  applies  for  the  case  L  ^  oo.  Slip  plane  6  now  consists  of  the 
two-sided  flat  region  {9  =  ±7r,  vq  <  r  <  Ro).  Body  05  is  multiply  connected  when  it  contains  the  slip  plane  (i.e.,  when 
it  is  treated  as  an  annular  tube)  since  9  is  multi-valued.  Domain  5B\©  is  simply  connected,  as  is  its  image  in  the  reference 
configuration  05o\6o. 

Equilibrium  equations  are  solved  on  05\6.  From  (3.19),  in  this  domain  the  strain  energy  function  in  (2.53)  becomes 

W  =  W{F)  =  W[D{F)]=W[E{F)].  (3.27) 

As  defined  in  Appendix  C  (See  Supplementary  Material,  Wiley  online  library),  Eulerian  and  Lagrangian  strain  measures 
(both  referred  to  material  coordinates  {AT'^})  are,  respectively, 

£)  =  i(l  -  F=i(FTF-l).  (3.28) 

Cauchy  stress  is,  from  (C.6), 

a  =  J-ipF^  =  J-^F^'^SF-^  =  J'^F^F^,  (3.29) 

where 

P  =  dWldF,  S  =  dW/dD,  S  =  dW/dE.  (3.30) 

Static  equilibrium  equations  are 

V-<T  =  0,  Vf,cr“'’  =  0.  (3.31) 

Internal  boundary  conditions  are 

|m]  —  b  on  6;  f|r  =  ro  =  =  per  a^’’\r  =  ro  =  —p-  (3.32) 

The  first  of  (3.32)  prescribes  a  displacement  jump  equal  to  the  true  Burgers  vector  across  the  slip  plane,  as  in  (3.9).  The 
second  of  (3.32)  assigns  a  traction  t  at  the  core  surface  corresponding  to  a  constant  radial  Cauchy  pressure  p,  with  a  unit 
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vector  in  the  radial  direction  as  defined  in  Appendix  A  (See  Supplementary  Material,  Wiley  online  library).  Traction  along 
the  core  must  be  a  uniform  pressure  to  maintain  consistency  with  the  symmetry  inherent  in  ansatz  (3.1).  External  boundary 
conditions  are 


r  =  R.o 


z=±Ll2 


=  0  or  s(z)  ^0  as  L 


(3.33) 


The  first  of  (3.33)  prescribes  traction  free  conditions  on  the  surface  of  the  cylinder  at  r  =  Rq.  The  second  of  (3.33) 
prescribes  either  null  average  axial  force  for  a  cylinder  of  finite  length,  or  no  average  axial  strain  for  an  infinite  cylinder. 

An  iterative  solution  procedure  in  Eulerian  coordinates  is  invoked,  similar  to  that  of  Teodosiu  and  Soos  [81,  82].  The 
primary  difference  between  the  present  isotropic  nonlinear  elastic  analysis  of  the  screw  dislocation  and  that  in  [81]  is 
that  here  a  different  (Eulerian)  strain  energy  potential  W  is  used,  in  contrast  to  the  Lagrangian  strain  energy  potential 
W  used  in  [81,  82].  A  second  difference  is  that  possible  axial  strains  are  considered  for  a  dislocation  of  finite  length,  in 
contrast  to  [81]  that  considered  only  the  infinitely  long  dislocation  and  imposed  s  =  0.  In  an  earlier  analysis  of  the  screw 
dislocation  in  isotropic  nonlinear  elasticity,  Seeger  and  Mann  [72]  used  an  iterative  procedure  in  Lagrangian  coordinates, 
with  a  Lagrangian  strain  energy  function,  and  allowed  for  axial  strain,  but  did  not  address  core  pressure.  The  general  method 
of  solution  in  Lagrangian  coordinates  was  described  formally  by  Gairola  [38]. 

According  to  the  method  of  solution  [81, 82],  spatial  displacement  u  is  written  in  the  series 

u{x)  =  kuo{x)  +  k^ui{x)  +  •  •  •  ,  (3.34) 


where  A:  >  0  is  a  small  scalar  parameter  whose  particular  value  will  not  affect  the  final  solution.  Internal  boundary  condi¬ 
tions  are  of  the  assumed  form 


b  =  kb,  p  =  kp,  (3.35) 

which  is  justified  by  the  vanishing  of  b  with  p  [81, 82].  Truncating  (3.34)  at  second  order"^  [i.e.,  omitting  terms  of  0(fc^)] 
and  prescribing  a  strain  energy  potential  of  either  of  the  latter  two  forms  in  (3.27)  results  in  a  series  expression  for  Cauchy 
stress: 


a  =  kao  +  k^ai+T),  ag"  =  (3.36) 

where  are  second-order  elastic  constants  (equivalent  to  those  introduced  in  Appendix  C  (See  Supplementary  Material, 
Wiley  online  library)),  and  where  t  depends  on  Vito  and  second-  and  third-order  elastic  constants.  Substituting  into  (3.31), 
two  sets  of  equilibrium  equations  result: 

V  •  CTo  =  0,  V  •  <Ti  =  —V  •  T.  (3.37) 

The  first  of  these  equations  is  solved  in  conjunction  with  internal  boundary  conditions  linear  in  k,  i.e.,  (3.35),  providing  a 
solution  for  Ug .  The  second  of  these  equations  is  then  solved  for  Ui ,  applying  continuity  of  Ug  across  6  and  null  traction  at 
r  =  rg,  along  with  the  previously  obtained  solution  for  Ug  to  be  used  to  compute  t.  Eor  both  solutions,  vanishing  traction 
at  r  =  i?o  is  imposed. 

The  general  method  of  solution  is  now  applied  to  the  screw  dislocation,  where  W{D)  is  used  for  the  strain  energy 
potential,  making  the  present  analysis  and  results  different  from  others  mentioned  already  [38,72,81].  Displacement  (3.34) 
applied  to  (3.8)  results  in 

u{r)  =  kug{r)  +  k‘^ui{r),  w{<p)  =  kwg{<p)  =  kbip/{2Tr),  s{z)  =  ksg{z)  +  k'^  Si{z),  (3.38) 

where  w  is  necessarily  truncated  at  first  order  to  satisfy  (3.7)  and  (3.35),  with 

b=kbGz,  b  =  kb.  (3.39) 

With  respect  to  curvilinear  cylindrical  bases  {  Ga  ,  9“  } ,  the  inverse  deformation  gradient  matrix  on  Q3\  6  is 


'r' 

0 

0  ’ 

1  —  kug  —  k’^u'i 

0 

0 

II 

7 

0 

0' 

0 

= 

0 

1 

0 

0 

deZ 

dzZ 

0 

—kb/{2'K)  1  —  ks'g  —  k"^  s[ 

4 


Incorporation  of  terms  cubic  in  k  would  result  in  a  problem  depending  on  fourth-order  elastic  constants,  which  cannot  be  measured  completely,  and 
are  usually  unknown  for  real  crystals. 
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As  remarked  already,  the  domain  of  analysis  excludes  the  core  and  slip  plane  and  thereby  avoids  singularities  and  consider¬ 
ation  of  .  In  physical  components,  i.e.,  orthonormal  bases  {sa,  Ea},  (A. 36)  becomes,  for  the  isolated  screw  dislocation. 


1  —  kuQ  —  k'^u'i 
0 
0 


0 


1  —  kuQ  jr  —  k^uijr 
—kb/{2'Kr) 


0 

0 

1  —  fcsg  —  k‘^s[ 


(3.41) 


noting  that  \fGjg  =  Rjr  =  \  —  ujr.  Henceforth,  such  physical  components  will  be  used  for  analysis  of  the  screw 
dislocation.  Volume  ratio  is 


J  ^  =  det[F  =  (1  —  kuQ  —  k^u'i){l  —  kuo/r  —  fc^ui/r)(l  —  /csq  —  k‘^s[) 

=  1  —  Uq  jv  4”  5q)  -f  0(A:^ ) . 


(3.42) 


The  next  step  is  computation  of  components  of  strain  tensor  D  in  terms  of  the  quantities  {k,  b,  u' ,  ujr^  s').  Subsequent  alge¬ 
bra  becomes  analytically  intractable  when  all  such  quantities  are  retained  in  the  general  solution.  Therefore,  two  particular 
cases  are  considered  in  what  follows: 


p  =  0=>so  =  0;  or  L— >oo=^s  =  0.  (3.43) 

The  first  case  is  consistent  with  a  priori  knowledge  that  in  the  linear  elastic  solution  of  a  screw  dislocation,  axial  strain 
vanishes  when  core  pressure  vanishes^.  This  is  the  Eulerian  analog  of  the  form  of  Lagrangian  displacement  function  s(Z) 
assumed  by  Seeger  and  Mann  [72].  The  second  case  in  (3.43)  is  consistent  with  the  second  set  of  boundary  conditions  in 
(3.33),  and  leads  to  the  condition  that  all  displacements  are  independent  of  z,  as  should  be  the  case  for  an  infinitely  long 
dislocation  line.  Therefore, 

s(z)  =  k^si(z)  (3.44) 

is  used,  with  si  =  0  invoked  when  boundary  conditions  corresponding  to  the  second  case  are  imposed. 

In  physical  cylindrical  components,  strain  tensor  D  is,  from  (3.28)  and  (3.41), 

D  —  DjiifEii  ®  Eji  +  DqqEq  0  Eq  +  DzzEz  ®  Ez  +  Dqz{Eq  ®  Ez  +  Ez  ®  Eq);  (3.45) 

Drr  =  kuQ  +  k'^{u[  -  Dqq  =  kuo/r  +  k'^{ui/r  -  k(uo/r)^], 

Dzz  =  k'^Wi  -  b'^ Dqz  =  Dze  =  kb/{4TTr)  -  k^buo/iAnr'^),  (3.46) 

where  terms  0{k^)  have  been  truncated.  Strain  energy  potential  per  unit  reference  volume  (C.7)  is  used,  truncated  at  third 
order: 

W  =  \CabcdDabDcd  +  \CabcdefDabDcdDef 

=  ^CapDaDjj  +  kCaPjDaDpD^,  (3.47) 

where  Greek  subscripts  denote  Voigt  notation.  For  an  isotropic  material,  second-order  constants  Cap  depend  on  (A,  /i)  from 
(C.ll),  and  third-order  constants  Cap-y  depend  on  (f'l ,  f’2  >  t'a )  from  (C.12).  Stress  S  —  dW  jdD  is,  omitting  terms  of 
0(9), 

S  =  SrrEr  (g)  Er  +  SqqEq  (g  Eq  +  SzzEz  (g  Ez  +  Sqz(Eq  (g  Ez  +  Ez  (g  Eq)\  (3.48) 


Srr  =  fc{(A  -I-  2/r)ito  -f  Xuq/t} 

+  fc^{(A  -I-  29)(u'i  —  Uo^/2)  -I-  A[Mi/r  —  |(ito/r)^  +  —  &^/(87r^r^)] 

-I-  4(i>i  -|-6i>2  +  Svi)u'^  +  (z>i  +  2v2)[\(uo/rY  +u'QUQ/r]  +  V2b’^ / (S9 9)) ,  (3.49) 

Sqq  =  fc{(A  -I-  29)uolr  +  Aug} 

-I-  9{{X  +  2fj,)[ui/r  -  |(■uo/r)^]  +  A[m(  -  Uq  /2  +  s(  -  b‘^/{899)] 

+  k(i)y  +  Qj>2  +  8f’3)('Uo/r)2  -p  (j>i  -p  2£'2)(moV2  -f  UgUoA) 

+  (E2  +  2B3)by(899)},  (3.50) 


5 


It  will  be  shown  explicitly  later  how  positive  core  pressure  p  is  linearly  proportional  to  axial  contraction  in  the  linear  solution. 


(c)  2013  WILEY- VCH  Verlag  GmbH  &  Co.  KGaA,  Weinheim 
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Szz  =  HKu'a  +  uo/r)} 

+  fc^{(A  +  2fi)[s[  -  6^/(87r^r^)]  +  A[ui  +  Ui/r  -  u'^ /2  -  |(ito/r)^] 

+  |(j>i  +  2v2)[uq  +  {uo/rY]  +  PiUgUo/r  +  (z>2  +  2v-i)b‘^  /  {8tt^  )} ,  (3.51) 

Sez  =  Sze  =  k{fib/{2TTr)}  -  /c^{^&uo/(27rr^)}.  (3.52) 

Using  (3.41)  and  (3.42),  Cauchy  stress  a  =  is,  omitting  terms  of  0{k^), 

a  =  (Trrer  0  6,.  +  aggee  ®eg+  +  ag^^^eg  0  0  eg)-,  (3.53) 

(Trr  =  fc{(A  +  2^)Uq  +  Alto/r} 

+  /c^{(A  +  2fi)u[  +  X{ui/r  +  s'l)  —  |(7A  +  14/r  —  Oi  —  6t>2  —  8z>3)uq^ 

-  (4A  +  2^  -  z>i  -  2i'2)uQUolr  -  i(3A  -  Oi  -  2z>2)(uoA)^ 

-  (A  -  i>2)&V(87tV^)},  (3.54) 

o-gfl  =  A:{(A  +  2^)uo/r  +  Amq} 

+  fc^{(A  +  2^)ui/r  +  A(u'i  +  s'l)  —  |(3A  —  Oi  —  2v2)u'i} 

-  (4A  +  2^  —  vi  —  2i>2)MoUo/r  —  |(7A  +  14^  —  vi  —  6z>2  —  8i>3)(uo/r)^ 

-  (A  +  2^  -  t>2  -  2t>3)6V(87rV2)},  (3.55) 

cTzz  =  fc{A(Mo  +UQ/r)} 

+  fc^{A(it'i  +  ui/r)  +  (A  +  2^)s'i  —  A(3A  —  z>i  —  2t>2)Mo^ 

-  (2A  -  z>i)uoMo/r  -  |(3A  -  Oi  -  2v2){ugi/rY 

-  (A  +  4/r  -  i>2  -  2t>3)6^/(87rV^)},  (3.56) 

ogz  =  o^g  =  fc{/i6/(27rr)}  -  k‘^{{X  +  fi)buQ / {2nr)  +  (A  +  4/r)5uo/(27rr^)}.  (3.57) 

These  are  consistent  with  (3.36).  From  symmetry  of  Cauchy  stress,  applying  (A.39)  in  spatial  components,  equilibrium 
equations  V  •  cr  =  0  reduce  to 

drarr  +  {a-rr  -  crgg)/r  =  0,  dgagg  f r  +  d^ag^  =  0,  dgag^ /r  +  =  0-  (3.58) 

First  consider  equilibrium  conditions  corresponding  to  stress  components  linear  in  k,  i.e.,  the  first  of  (3.37),  where 

<’’o(?')  =  [(A  +  2fi)uQ  +  Xuo/rje,.  0  6^  +  [(A  +  2^)uo/r  +  Xu'Q]eg  (g)  eg 

+  [A(uo  +  UQ/r)]e^  ®  +  [^ib / {2'Kr)]{eg  ®  e^  +  ®  eg).  (3.59) 

Since  (Tq  does  not  depend  on  9  or  z,  the  second  and  third  equilibrium  equations  in  (3.58)  are  trivially  satisfied  to  first  order 
in  k.  The  first  order  radial  equilibrium  equation  becomes 

(A  +  2^){uq  +  Xu'q/t  —  uo/r^)  =  0.  (3.60) 

This  is  a  homogeneous  second-order  ordinary  differential  equation  of  Cauchy-Euler  type,  with  general  solution 

Mo  =  cir  +  C2/r.  (3.61) 

Applying  boundary  conditions  (cro)rr  |r=ro  =  ~P  ™d  (CTo)rr  |r=ito  =  the  solution  is 

Cl  =  x(l  -  2v),  C2  =  x^o,  X  =  b/(2l^)][fo/(^o  -  ^’o)]:  (3-62) 
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with  Poisson’s  ratio  v  defined  in  (C.14).  With  uq  now  known,  all  components  of  <To  can  be  found  using  (3.59).  To  first 
order  in  k,  radial  displacement,  radial  stress,  circumferential  stress,  and  axial  stress  vanish  if  core  pressure  is  omitted 
(X  =  0)^  whereby  the  solution  degenerates  to  the  usual  anti-plane  shear  description  of  a  screw  dislocation  in  isotropic 
linear  elasticity.  It  was  assumed  in  (3.44)  that  so('2^)  =  0.  If  this  restriction  is  relaxed  in  the  linear  solution,  then  the  axial 
equilibrium  equation  yields  (A  +  2p.)sQ  =  0  =>  sq  =  c^z.  If  the  average  axial  stress  (<70)2^  is  set  to  zero  over  any  cross 
section  of  the  cylinder  at  constant  z,  then  sq  =  —  2^) /(I  —  v)]xz  is  the  axial  contraction  that  occurs  in  conjunction 

with  positive  core  pressure  for  ^  >  0.  Since  Sq  =  0  when  x  =  the  first  assumption  in  (3.43)  is  now  verified. 

Now  consider  equilibrium  equations  for  stress  components  quadratic  in  k,  i.e.,  the  second  of  (3.37).  From  (3.54)-(3.57), 

<Ti(r,  z)  =  [(A  -f  2^)u[  +  X{ui/r  +  sj)]er  0  -f  [(A  -f  2fj,)ui/r  +  X{u[  +  s[)]ee  (g)  eg 

+  [(A  -I-  2/r)s'i  -I-  X{u\  +  Mi/r)]ez  ®  e^,  (3.63) 

T(r)  =  {x^(l  -  2j2)^(-9A  -  9/r  -I-  2vi  +  6£'2  -f  dDs) 

+  -  2:/)(4A  -f  14p  -  4j>2  -  8i>3)  -  (A  -  V2)b‘^ / 

+  -  5p  -I-  2v2  +  4i>3)]}er  ® 

+  {x^(l  -  2vy{-9X  -9^1  +  2i>i  -I-  6£'2  -I-  4j>3) 

+  il/r^)[x^Rl{l  -  2iy){-AX  -  14p  +  4z>2  +  8ps) 

-  {X  +  2fj,  -  i>2  -  2i>3)&^/(87r^)] 

+  (lA^)[x^-Ro(-^  -  5p  -f  2i>2  -f  ^h)]}eg  (g)  eg 

T  {x^(i  ~  2i^)^(— 5A  -f  2vi  -f  21)2)  -f  (l/r^)[(— A  —  4p.  V2  2v-,{)b’^ I (Stt^)] 

+  (lA^)[x^-Ro(-^  +  2i>2)]}e2  (g) 

-  {(1/0 [(1  -  2^)(2A-f  5/r)6x/(27r)] 

-f  {l/r^)[8^JibxRl/{2.^T)]}{eg  ®  e^  +  e^  ®  eg).  (3.64) 

The  radial  component  of  second-order  equilibrium  equation  in  (3.37)  is 

dr  (cTi  -f  T)rr  +  [(cTl  +  t)„.  -  (cti  -f  T)gg]/r  =  0.  (3.65) 

Substituting  from  (3.63)  and  (3.64),  this  becomes  an  inhomogeneous  second-order  ordinary  differential  equation  of  Cauchy- 
Euler  type: 


(A  -f  2fj,){u”  +  u[/r  -  Mi/0)  =  (l/0)(i>2  -f  i>3  -  A  -  /r)O/(47r0 

+  (1/0)(2z>2  +  4:>3  -  a  -  5p)4x'i?0 


(3.66) 


The  homogeneous  general  solution  to  this  equation  is  of  the  form  mih  =  Cir  -f  C'2/r.  The  particular  solution  uip  can  be 
found  using  Lagrange’s  method  (i.e.,  variation  of  parameters).  The  total  general  solution  is  the  sum  ui  =  mih  +  wip: 


,  C'z  ,  &^(A-f /r- i>2  -  i>3)  Inr  x^ (2i>2  +  4f'3  -  A  -  5^) 
«,(r)  =  C.r  +  -  +  — +  - T  + - 2(aTW - 


(3.67) 


Consistent  with  previous  arguments,  applying  vanishing  second-order  traction  end  conditions  (cri)rr-  =  —Trr  at  r  =  ro  and 
r  =  i?()  gives 


/r(A  +  H-P2  -P3)  ln{Ro/ro) 


Stt^  [  {X  +  fj,){X  +  2fi)  Rl  —  rl  j 


—  Z/Si 


xV 


-  [(1  —  2v){9X  +  9fi  —  2Pi  —  6z>2  —  4F3) 


2(A-f  p)2 
+  2(A  +  5p-2:>2-4i)3)(i?^/rg)], 


(3.68) 
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Co  = 


Stt^ 


i>2  -  A 

2fi 


+  (A  +  /r  -  j>2  -  O3)  —  + 


1  rg  In  i?o  —  i?o  In  rg 


2/r  (A  +  2^)(i?2-rg2);j 


(1  -  2r2)(4A  +  14^  -  4z>2  -  8O3)  (A  +  5^  -  21>2  -  4j>3)(i?o  +  rg) 


(3.69) 


2/i  2(A  +  2/r)rQ 

Since  (Ti  and  r  do  not  depend  on  6  and  since  {(Ji}ez  =  Tez  =  0,  the  second-order  circumferential  equilibrium  equation  in 
(3.58)  is  trivially  satisfied.  The  second-order  (i.e.,  proportional  to  fc^)  axial  equilibrium  equation  is 


^^[(cri)^^  +  Tzz]  =  0  ^  (A  +  2/r)s"  =  0. 

Its  general  solution  is,  upon  imposing  si  (0)  =  0  to  fix  translational  invariance, 
Si  =  C3Z. 


(3.70) 


(3.71) 


Therefore,  s[  =  C3  =  constant,  and  Cauchy  stress  is  of  the  form  a  =  cr(r).  For  an  infinite  dislocation  line  (L  00), 
the  second  of  (3.43)  requires  C3  =  0.  For  the  first  condition  in  (3.43),  corresponding  to  a  dislocation  of  finite  length  with 
p  =  0,  applying  the  null  axial  force  condition  over  any  cross  section  normal  to  the  dislocation  line: 


/TT  pR(i 

/  [(^1). 

-TT  J rn 


zjrdrdff  =  0 


'  —IT  Jro 

results  in  an  algebraic  equation  for  C3  that  yields 


C3  = 


62 

Tn(i?g/rg)' 

A  +  // 

47r2 

[  \ 

^(3A  -l-  2/i) 

A(i>2  +  O3  -  A  -  p)  ,  .  ,  .  .  A/r(A  +  p-£'2 -O3) 

- h  A  +  4/r  —  V2  —  21/3  — 


A  4“  2/i 


( A  +  /i)  (A  +  2/i) 


(3.72) 


(3.73) 


Recall  that  (cro)z2  =  0  when  p  =  0,  so  that  the  total  axial  force  from  terms  proportional  to  k  and  vanishes  in  this  case  as 
well.  Since  =  C^,  solutions  for  Ci  and  C2  are  complete  and  ui  (3.67)  is  fully  known,  as  is  second-order  Cauchy  stress 
<Ti  in  (3.63). 

Combining  first-  and  second-order  solutions  for  displacement,  (3.38)  becomes 


u{r)  =  k 


X(1  -  2i/)r  + 


xRl 


Cir  H - 


6^ (A  +  /i  -  O2  -  i>3)  In r  x^Roi‘^^2  +  4z>3  -  A  -  5/i) 
87r2(A  +  2/i)  r  2(A  +  2/i)r3 


(3.74) 


w{6)  =  kb9/{2Tr),  9  €  (— 7r,7r); 

s(z)  =  k^C^z. 


(3.75) 

(3.76) 


Recall  from  (3.35)  that  kb  =  b  and  kp  =  p,  and  from  (3.62)  that  P-  terms  in  the  solution  linear  in  k  are  proportional 
to  b  or  X-  All  terms  in  the  solution  quadratic  in  k  are  proportional  to  b^  or  x^  •  Therefore,  k  does  not  affect  the  final  solution. 
Combining  first-  and  second-order  solutions  for  Cauchy  stress,  (3.36)  becomes 


a{r)  =  kcTo  (r)  +  k^  [cti  (r)  +  T(r)] , 


(3.77) 


where  (Tg  is  given  by  (3.59)  with  all  terms  linear  in  b  or  x,  and  cti  and  t  are  given  by  (3.63)  and  (3.64),  respectively,  with 
all  terms  quadratic  in  6  or  x  or  bilinear  in  bx- 

Strain  energy  per  unit  reference  volume  W  of  (3.47)  can  be  written 

W  =  iPag  :  Do  +  k^ao  :  Do  abDo  cdDo  ef  +  0{k^),  (3.78) 

where  Dg  is  the  first-order  (linear  in  k)  part  of  strain  in  (3.45): 

Do  =  u'oEji  Cl  Ee,  +  {uo/r)EQ  C  Eq  +  [&/(47rr)](£le  C  Ez  +  Ez  C  Eq),  .  (3.79) 
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and  k^Di  =  D-  kDo  is  the  second-order  part  of  strain.  Considered  here  is  the  infinitely  extended  dislocation  line.  The 
first  term  in  (3.78)  is  computed  as 

|tTo  :  Do  =  fib^ +2/xx^[(l  -  2iy)  +  Ro/r'^].  (3.80) 


Assuming  contributions  from  higher-order  products  of  radial  displacement  and  its  derivatives  are  small  compared  to  those 
involving  shear,  second  and  third  terms  are 


o-Q  :  A  ~ 


1  /^ABCDEF  7~)  7~)  7~) 

wLy  A B  ^0  CD  -^0  EF 


62 


X 


(1  -  2v){v2  +  h)  + 


RC 


(3.81) 

(3.82) 


Strain  energy  per  unit  current  volume  kW  /  J  where  J  ^  is  given  by  (3.42);  integrating  this  over  the  cross  section  gives 
the  energy  per  unit  length  T',  where  x  =  ^X- 


p0  =  -K  nr  =  B.Q 

T'  =  /  /  J~^Wrdrd9 

J 6  =  —'k  J r=ro 

«  A  —  {m[1  -  2x(l  -  ‘2i')]  +2x(l  -  2i/)(z>2  +  £'3-  X-  fi)} 

ATT  ro 

+  _  2^)2  ^  _  2^(1  _  2^)] 

+  2Trx{Rl/rl  -  l){b\^3  -  £t)/i8^)  +  xkiRl[l  -  2x(l  -  2iy)]}.  (3.83) 


The  linear  elastic  energy  per  unit  length  is  recovered  by  setting  b^XtX^  0^  when  core  pressure  vanishes,  this  linear 
solution  is  the  conventional  T'  =  / {4:7:)]  ln(i?o/ro). 


3.3  Solution  for  ideal  crystal  with  Cauchy  symmetry 

The  complete  solution  depends  on  9  parameters:  geometric  constants  {b,ro,Ro),  elastic  constants  {X,  i/2,1'3),  and 

core  boundary  condition  p.  The  solution  is  further  examined  for  a  reduced  set  of  parameters.  Choosing  6  =  Tq  [56],  the 
geometry  and  core  pressure  become  respective  functions  of  dimensionless  parameters  A  and  po : 

A  =  Ro/ro  =  Ro/b>  I,  po=p/K.  (3.84) 

Imposing  Cauchy  symmetry  on  the  elastic  constants,  which  is  a  physically  reasonable  assumption  for  alkali  metals  and  ionic 
crystals  such  as  the  alkali  halides  [40],  produces  relations  (C.20)  and  (C.22)  of  Appendix  C  (See  Supplementary  Material, 
Wiley  online  library).  The  number  of  elastic  constants  is  now  two,  the  linear  bulk  modulus  K  and  pressure  derivative  of 
tangent  bulk  modulus  K'  =  Bq'. 

X  =  p=lK,  =  ^>2  =  -f  (|A-6),  i>3  =  -f(fA-9).  (3.85) 

Thus  1/  =  \  and  all  other  elastic  constants  are  proportional  to  K.  Typical  values  for  K'  in  natural  crystalline  solids  are 
2  <  K'  <  7  [45],  with  AT'  «  4  a  common  approximation  if  high  pressure  data  are  unavailable.  The  reduced  set  of 
5  parameters  entering  the  solution  is  now  {b,A,po,K,K').  Noting  that  x  =  \po/{A?  —  1)  and  that  all  coefficients  of 
displacement  depend  only  on  dimensionless  ratios  of  elastic  constants,  normalized  displacement  u/b  depends  only  the  set 
of  3  parameters  {A,po,  K').  Since  stress  components  and  strain  energy  are  linearly  proportional  to  an  elastic  constant, 
dimensionless  stress  and  energy  a/K  and  4/  jK  also  depend  only  on  this  set.  It  was  verified  fhat  fhe  dimensionless  solution 
does  nol  depend  on  b  or  its  units. 

Normalized  radial  displacement  u/rg  =  u/b  is  shown  versus  Eulerian  radial  coordinate  r  in  Fig.  1.  In  this  and  all 
later  plots  of  results  versus  radial  coordinate,  the  abscissa  begins  at  its  minimum  valid  value  r/rg  =  1.  The  effect  of  core 
pressure  is  shown  in  Fig.  1(a)  for  a  fixed  cylinder  size  of  A  =  i?g/rg  =  100:  dilatation  (i.e.,  radial  expansion)  increases 
with  increasing  po  but  decreases  rapidly  with  increasing  r.  The  effect  of  nonlinear  elastic  constant  K'  is  shown  in  Fig.  1(b) 
for  a  fixed  core  pressure  po  =  0.1  and  and  an  infinitely  long  dislocation,  i.e.,  fixed  ends  {C3  —  0):  dilafation  increases 
with  increasing  K' .  The  effect  of  K'  is  shown  in  Fig.  1(c)  for  null  core  pressure  but  average  free  axial  conditions  (3.72): 
dilatation  increases  in  magnitude  with  increasing  K' ,  is  positive  for  K'  >  4,  and  is  negative  for  K'  <  4.  Values  of  C3 
computed  from  (3.73)  are  shown  in  the  legend;  recall  that  C3  >  0  corresponds  to  axial  extension,  while  A  <  0  corresponds 
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5  10  15  20  25  30 


r/r^ 


(a) 


(b) 


r/to 


(c) 


(d) 


Fig.  1  Total  radial  displacement  for  screw  dislocation:  (a)  variable  core  pressure  po  (b)  variable  K' ,  po  =0.1  (c)  variable 
K' ,  Pa  =  0  (d)  variable  cylinder  size  Rq  =  Aro . 


to  axial  contraction.  For  K'  =  4,  u  ^  0  for  r  >  vq,  suggesting  that  core  pressure  rather  than  elastic  nonlinearity  may  be  the 
primary  contributor  to  dilatation  in  a  typical  crystalline  solid.  Recall  in  the  linear  elastic  solution  po  =  0  =>  mq  =0.  The 
effect  of  cylinder  size  Rq  is  shown  in  Fig.  1(d):  for  A  >  100,  radial  displacements  are  indistinguishable,  while  for  A  <  10, 
increased  dilatation  is  observed  since  the  free  surface  is  closer  to  the  core. 

Normalized  Cauchy  stress  field  a  jK  is  shown  in  Fig.  2  for  A  =  100  and  K'  =  4.  Radial  stress  is  shown  in  Fig.  2(a); 
the  vertical  axis  is  truncated  at  (Jrr  =  — O.SAT  since  arr  =  —Kpo  at  r  =  rg.  Maximum  magnitudes  at  r  =  Tq  are  on  the 
order  of  |po  I  -  Radial  and  circumferential  stress  components  are  of  opposite  sign  and  decay  rapidly  within  r  <  15ro.  Axial 
components  are  negative  close  to  the  core  and  slightly  positive  farther  away  from  the  core,  and  decrease  rapidly  within 
r  <  5ro.  All  stresses  converge  towards  a  single  curve  for  each  component  for  r  >  20ro  that  corresponds  to  the  linear 
elastic  solution  with  vanishing  core  pressure.  For  large  r,  normal  stresses  decay  to  zero  and  shear  stress  ag^  [Fig-  2(d)]  is 
proportional  to  1/r. 

Nonlinear  and  linear  solutions  for  stress  and  energy  per  unit  length  are  compared  in  Fig.  3  and  Table  1,  respectively. 
Stress  components  shown  in  Fig.  3  correspond  to  fixed  parameters  A  =  100, po  =  0.1  (linear  and  nonlinear)  and  K'  =  4 
(nonlinear).  Radial  and  circumferential  stresses  obtained  from  the  nonlinear  solution  are  nearly  indistinguishable  from 
corresponding  stresses  from  the  linear  solution.  Axial  stress  is  negative  and  larger  in  magnitude  for  the  nonlinear  solution 
compared  to  the  linear  solution  for  r  <  drg.  Shear  stress  is  of  the  same  positive  sign  but  smaller  in  magnitude  for  the 


Table  1  Strain  energy  /{tvRIK)  x  10^  for  screw  dislocation,  A  =  100. 


PO 

o 

II 

K'  =  2 

K'  =  4 

K'  =  6 

K'  =  8 

linear 

1.0 

0.929 

0.918 

0.906 

0.896 

0.885 

0.903 

0.1 

0.809 

0.798 

0.787 

0.776 

0.765 

0.783 

0.0 

0.700 

0.700 

0.700 

0.700 

0.700 

0.700 

-0.1 

0.758 

0.769 

0.780 

0.791 

0.801 

0.783 
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5  10  15  20  25  30 


r/rg 


(a) 


(b) 


(c)  (d) 


Fig.  2  Cauchy  stress  field  for  screw  dislocation  with  /f'  =  4  and  A  =  100:  (a)  radial  stress  (b)  circumferential  stress  (c) 
axial  stress  (d)  shear  stress. 


Fig.  3  Comparison  of  nonlinear  and  linear  elastic  solutions  for  stress  field  of  screw  dislocation  with  A  =  100,  po  =0.1 
and  K'  =  4. 


nonlinear  solution  for  r  <  4ro .  The  total  stress  field  from  the  nonlinear  solution  becomes  indistinguishable  from  the  linear 
solution  for  r  >  lOrg.  Solutions  for  normalized  strain  energy  in  Table  1  are  obtained  using  (3.83),  where  terms  oc  are 
omitted  in  the  linear  solution.  For  pg  =  0,  K'  does  not  affect  the  solution,  and  nonlinear  and  linear  predictions  coincide. 
For  other  values  of  core  pressure,  energy  predicted  using  the  nonlinear  solution  may  exceed  the  linear  solution  depending 
on  the  value  of  K' . 
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4  Wedge  disclination 

4.1  Geometric  description 

Consider  a  deformed  body  05  with  boundary  i9®.  In  the  reference  configuration,  the  image  of  this  body  and  its  boundary  are 
05o  and  5)Bo.  Notation  for  cylindrical  and  Cartesian  coordinate  systems  is  the  same  as  that  used  in  Sect.  3  and  Appendix  A 
(See  Supplementary  Material,  Wiley  online  library). 

An  Eulerian  description  of  deformation  is  invoked,  similar  to  that  of  Vladimirov  et  al.  [86],  though  a  different  strain 
energy  potential  and  core  boundary  conditions  will  be  used  here.  Let  the  unit  tangent  line  of  the  disclination  be  oriented 
along  the  z-axis,  i.e.,  ^  =  g^-  The  general  ansatz  for  deformation  due  to  a  wedge  disclination  in  an  isotropic  body  is 

R{r)=r  —  u{r),  Q{9)  =  9  —  ilnp/{2'K),  Z{z)  =  z  —  s{z).  (4.1) 

Radial  displacement  is  u{r),  and  a)  is  a  constant  associated  with  the  Frank  vector.  When  w  >  0,  the  wedge  disclination  is 
labeled  “positive”,  corresponding  to  removal  of  a  wedge  of  material  and  tensile  circumferential  stress.  When  w  <  0,  the 
disclination  is  “negative”,  corresponding  to  insertion  of  a  wedge  of  material  and  compressive  circumferential  stress.  Angle 
9  is  multi-valued  upon  complete  revolutions  around  the  disclination  line,  and  ip  G  (— tt,  tt)  is  defined  in  (3.2).  Surface  6, 
where  0(9)  is  discontinuous,  corresponds  to  the  half-plane  (y  =  0,x  <  0),  or  9  =  zLtt.  On  1B\©,  the  second  of  (4.1)  can 
be  inverted  to  9  =  k0,  where  k  =  2Trl(2'K  —  uj). 

Displacement  u  =  x  —  JC  in  (A. 38)  in  physical  curvilinear  basis  {ea(9)}  is 

u  =  {r  —  (r  —  u)  cos[a)(p/(27r)]}e,.  +  {(r  —  u)  sin[a)(p/(27r)]}eg  +  se^.  (4.2) 

Displacement  is  continuous  except  for  that  due  to  the  jump  in  p  across  6: 

|01  =  0|^=, -0|e=_,  =  -a;,  (4.3) 

M  =  M|e=7r  -  =  (u-  r) { |cos[(U(^/(27r)]] 6,.  -  |sin[wv?/(27r)]]e0}.  (4.4) 

For  small  u  and  small  tu,  this  reduces  to  the  usual  approximation  |i6]  «  ruies  [14,32]. 

In  curvilinear  coordinates,  the  inverse  deformation  gradient  from  (4.1)  is 

F-i  =  F-\^Gr  0  9’-  +  F-i®Ge  0  g'  +  Gz  0 

=  {dR/dr)Gji  0  g’’  +  [1  -  ii!dgp/{2TT)]GQ  0  g®  +  (dZ/dz)Gz  0  g"" 


=  [1  -  u'{r)]GR  0g’'  +  [1  -  u;d0p/{2TT)]Ge  0  g®  +  [1  -  s'{z)]Gz  0  gT  (4.5) 

Applying  the  generalized  derivative  from  (3.5), 

=  1  —  [Cj/(2-K)][—rs\n9dxP  +  rcos9dyp\  =  1  —  Cj/(2t:)  +  0rH(— x)(5(y).  (4.6) 

Define  fhe  funcfion  w(9)  and  ifs  derivative 

w  =  d>0/(27r),  w' =  dw/d0  =  a)/(27r).  (4.7) 

To  fix  translafional  invariance  of  the  origin  and  ensure  volume  remains  positive, 

u(0)  =  0,  s(0)  =  0;  (l-M/r)(l-u')(l-w')(l-s')  >  0-  (4-8) 

Physically,  constraint  |a)|  <  27r  also  applies. 

The  multiplicative  decomposition  of  lattice  deformation  in  the  first  of  (2.46)  is  now  applied: 

F  =  F^F\  p-t  ^  pt-i  pn-t  ^  (4  9) 

The  “elastic”  part  of  F^'^  is  defined  as  that  which,  in  curvilinear  coordinates,  is  continuous  and  single-valued  over  05: 

=  [1  -  u'(r)]GR  0  g’-  +  [1  -  zn'jGe  0  g®  +  [1  -  s\z)]Gz  0  gT  (4.10) 

It  follows  from  (4.5),  (4.6),  and  (4.9)  that 

=  LorR(-x)S(y)Ge  0  g®  +  (4.11) 
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Decomposition  (4.9)  implies  existence  of  a  locally  unstressed  intermediate  configuration  IB  for  each  material  element. 
Selecting  coincident  cylindrical  coordinate  systems  on  IB  and  *Bo,  such  that  ga{X)  =  S^Ga{X)  [15],  and  solving  for 

pi— 10  _  p— 10  pE0 

^  .0  —  ^  .$  ^.0  ’ 

=  {1  +  [27rrw/(27r  -  w)]H(-a:)(5(2/)}G0  (g)  G®.  (4.12) 

From  the  above  description, 

F^^{x)  =  F^-^  Va:G»\6,  F-\x)  =  F^-^  F^-'^  on  6.  (4.13) 

These  can  be  inverted  for  F{x).  Since  0(0)  is  multi-valued  across  ©,  the  referential  image  of  6  is  not  a  single  half-plane 
as  was  the  case  for  the  screw  dislocation  in  Sect.  3. 

The  disclination  density  is  now  considered.  A  spatial  density  tensor  6  and  Frank  vector  u)  of  the  following  forms  are 
consistent  with  geometry  of  the  problem: 

B  =  0^^gz®gz=Cu5{r)g:,®g^=5{r)u®^-,  u  =  ujg^.  (4.14) 

To  complete  the  description  of  kinematics,  variable  Q  entering  the  connection  (2.19)  is  sought  that  yields  (4.14)  when  used 
in  (2.28)  and  (2.30).  From  the  second  of  (2.33),  the  only  nonvanishing  covariant  components  of  curvature  in  are,  in  spatial 
Cartesian  coordinates, 

Rxyxy  —  Ryxyx  —  Rxyyx  —  ^yxxy  —  0  ■  (4.15) 


From  (2.28), 


Rxyxy  —  [xQy][xy]  +  Txy  Qa\xy])  Qx[xQ\y\y]a  +  Qy[xQ\x\y]aj 


(4.16) 


with  summation  on  index  a  and  (-([sj/]  =  |  [(Oxi/  ~  (Oi/x]-  Existence  of  a  field  Q  that  simultaneously  obeys  (4.14)-(4.16) 
and  all  of  (2.47)  is  unproven  at  present.  However,  the  following  approximation®  is  known: 


Q  =  (2jH{—x)6{y){ey  ®  By  ®  —  By  ®  ®  By),  (4.17) 

Qyyx  —  Qyxy  —  Ci.^H(  x)5{jj).  (4.18) 

Then  (4.16)  becomes 

Rxyxy  —  ^xQyxy  QyxQyxy  ~  ^xQyxy  —  x)S(y)  —  o)S(^t)  —  0  ,  (4.19) 

where  terms  of  0(a)^)  have  been  omitted  in  the  approximation,  yielding  a  result  consistent  with  (4.14)  and  (4.15).  In  the 
linear  approximation,  the  contribution  to  the  local  Burgers  vector  from  the  disclination  is,  with  n  —  g^  in  (2.48): 

db  =  :  {On  (g)  x)da  «  a)  x  a;da  =  rCjS{r)drdOgg .  (4.20) 


Integral  f  db  then  vanishes  identically  for  domain  a  with  constant  normal  g^  parallel  to  ^  when  the  disclination  passes 
through  r  =  0.  But  if  the  coordinate  system  is  translated  such  that  ^  intersects  planar  area  a  at  coordinates  {xq  = 
To  cos  00,2/0  =  ’’0  sin  0o),  then  the  disclination  density  becomes  0^^  =  u)5{x  —  X[))5{y  —  yg)  and  b  =  f  db  =  roUJBg(0o)- 


4.2  Nonlinear  elastic  analysis  and  general  solution 

A  nonlinear  elastic  boundary  value  problem  for  the  wedge  disclination  is  constructed  as  follows.  Let  body  05  consist  of  a 
cylinder  of  outer  radius  Rq  with  the  disclination  line  located  along  r  =  0.  A  cylindrical  core  region  r  <  rg  is  removed 
from  the  body  and  excluded  from  the  solution  domain.  The  length  of  the  cylinder  is  L,  such  that  the  body  is  contained 
within  the  limits  though  the  forthcoming  analysis  also  applies  for  the  case  L  ^  oo.  Plane  0  now  consists 

of  the  two-sided  flat  region  {9  =  =b7r,ro  <  r  <  Rq).  Body  ^  is  multiply  connected  when  it  contains  this  plane  (i.e., 
when  it  is  treated  as  an  annular  tube),  since  6  is  multi-valued.  Body  ^\0  is  simply  connected.  Its  image  in  the  reference 
configuration,  !Bo\0o»  is  simply  connected  when  a;  >  0,  but  contains  an  overlap  of  material  across  ©o  when  c2;  <  0. 

^  This  approximation  is  not  fully  consistent  with  the  second  of  (2.47)  since  T^'y  —  —  ^  Qy^  does  not  vanish,  leading  to  singular  dislocation  density 
a'y^  —  Qyx  on  ©.  Similar  aitifacts  exist  in  linear  representations  of  discrete  disclination  lines  [32].  The  result  is  inconsequential  here  since  the 
elastic  solution,  which  does  not  depend  explicitly  on  the  choice  of  Q,  is  sought  only  over  ^\©. 
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The  solution  procedure  parallels  that  of  the  screw  dislocation  of  Sect.  3.2.  Equilibrium  equations  are  solved  on  iB\6. 
From  (4.13),  in  this  domain  the  strain  energy  function  in  (2.53)  becomes 

W  =  W{F)  =  W[D{F)\=W[E{F)],  (4.21) 

where  W  is  used  in  the  solution  that  follows.  Equations  (3.28)-(3.31)  apply. 

Internal  boundary  conditions  are 


|0]  =  —Cj  on  6;  f|r  =  ro  =  =  per  tT’'’’|r  =  ro  =  —p- 


(4.22) 


The  first  of  (4.22)  prescribes  a  coordinate  jump  across  the  plane  of  discontinuity  as  given  in  (4.3).  The  second  of  (4.22) 
assigns  traction  t  at  the  core  surface  corresponding  to  constant  radial  Cauchy  pressure  p,  with  a  unit  vector  in  the 
radial  direction.  Traction  along  the  core  must  be  a  uniform  pressure  to  maintain  consistency  with  symmetry  in  ansatz  (4.1). 
External  boundary  conditions  are 


t^rdrdO  j 

ro  J 


z=±LI2 


=  0  or  s(z)  ^0  as  L 


(4.23) 


The  first  of  (4.23)  prescribes  traction  free  conditions  on  the  surface  of  the  cylinder  at  r  =  Rq.  The  second  of  (4.23) 
prescribes  either  null  average  axial  force  for  a  cylinder  of  finite  length,  or  no  average  axial  strain  for  an  inhnite  cylinder. 

An  iterative  solution  procedure  in  Eulerian  coordinates  is  invoked,  similar  to  that  of  Vladimirov  et  al.  [86].  The  primary 
difference  between  the  present  isotropic  nonlinear  elastic  analysis  of  the  wedge  disclination  and  that  in  [86]  is  that  here 
a  different  strain  energy  potential  W  is  used.  A  second  difference  is  that  possible  core  pressure  is  considered,  in  contrast 
to  [86]  that  treated  the  core  in  a  different  way.  Core  pressure  was  omitted  in  prior  linear  solutions  [30,32].  Equations  (3.34), 
(3.36),  and  (3.37)  apply.  Internal  boundary  conditions  are  of  the  assumed  form 


cP  =  ku>,  p  =  kp, 


which  is  justified  by  the  vanishing  of  u)  with  p.  Displacement  (3.34)  applied  to  (4.1)  and  (4.7)  results  in 
u{r)  =  kuo{r)  +  k‘^ui{r),  w{9)  =  kwo{6)  =  ku!6/{2Tr),  s{z)  =  kso{z)  +  fc^Si(z), 

where  w  is  necessarily  truncated  at  first  order  to  satisfy  the  first  of  (4.22),  with 


UJ  =  kujQz,  Cj  =  koj. 


With  respect  to  curvilinear  cylindrical  bases  {  G a  ,  } ,  the  inverse  deformation  gradient  matrix  on  03  \  6  is 


0 

1 

o 

1  —  fcrtg  —  k'^u'i 

1 - 

o 

o 

0 

0' 

0 

= 

0 

1  —  kwQ  0 

1 - 

o 

0 

Z' 

0 

0  1  —  fcsg  —  k‘^s[ 

In  physical  components,  i.e.,  orthonormal  bases  {sa,  Ea},  (A. 36)  becomes 

1  —  fcrtg  —  k'^u'i  0  0 

0  {1  —  kuQ  /  r  —  k"^  ui  /  r){l  —  kw'o)  0 

0  0  1  —  fcsg  — 

noting  that  \f(T[g  =  Rjr  =  1  —  ujr.  Volume  ratio  is 

J~^  =  (1  —  ku[^  —  k^u'i){l  —  ku^/r  —  fc^Mi/r)(l  —  fcwo)(l  ~  ~ 

=  1  —  fc(Wg  +  Mo/r  +  Wg  +  Sg)  +  0{k'^) . 

In  physical  cylindrical  components,  strain  tensor  D  is,  from  (3.28)  and  (4.28), 

D  =  Diij^Eii  (g)  E]i  +  DqqEq  (g)  Eq  +  DzzEz  ®  Ez] 


[F- 


{Aq 


(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 
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Drb.  =  ku'n  +  k‘^{u[  -  Dzz  =  ks^  +  k‘^{s[  - 

Dqq  =  k{wQ  +  uo/r)  +  k^[uilr  -  \{uQ/rY  -  \w'^  -  2w[)Uo/r],  (4.31) 

where  terms  0{k^)  have  been  truncated.  Strain  energy  potential  per  unit  reference  volume  (C.7)  is  used,  truncated  at  third 
order,  leading  to  (3.47)  which  also  applies  here.  Stress  S  =  dW/dD  is,  omitting  terms  of  O(fc^), 

S  =  SrrEr  (g)  Er  +  SqqEq  (g)  Eq  +  SzzEz  0  Ez;  (4.32) 

Srr  =  k{{X  +  2/i)ito  +  X{uo/r  +  Wq  +  4)} 

+  k^{{X  +  2fi){u[  -  Uq  /2)  +  X{ui/r  +  s[  -  \  [{uQ/rf'  +  w'^  +  -2'u;^uo/r) 

+  4(t>i  +  6t>2  +  +  4(t>i  +  2z>2)[(w^  +  uo/rf  +  s'i} 

+  2u'f^{uQ  /  r  +  w'q  +  So)]  +  i>iSo(uo/r  +  Wg)},  (4.33) 

Sqq  =  k{{X  +  2fi){wQ  +  uo/r)  +  X{uq  +  Sq)} 

+  A:^{(A  +  2fi){ui/r  -  |[(ito/?’)^  +  -  2wgUQ/r)  +  X{u[  -  v!^ /2  +  -  s'^  12) 

+  i(i>i  +  6t>2  +  8t>3)(w^  + 

+  \(y\  +  2v2)[u'q  +  Sq^  +  2(Mg  +  SQ)(mQ  +  uo/r)]  +  DiMqSq},  (4.34) 

Szz  =  k{{X  +  2/r)so  +  A(uo  +  Wg  +  Wg/r)} 

+  {(A  +  2fj){s[  -  S0V2)  +  A(ug  +  ui/r  -  4[ug2  +  (ug/r)^  +  w'^]  -  2wQUo/r) 

+  +  6f'2  +  8C'-^)s'q 

+  A(t>i  +  2i>2)['Uo^  +  {wq  +  uo/rf  +  2(wq  +  Mg  +  wgA)so]  +  t>iUo(M;o  +  wg/r)}.  (4.35) 

Using  (4.28)  and  (4.29),  Cauchy  stress  a  =  is,  omitting  terms  of  0{k^), 

a  =  fTrr-e,.  (g  e,.  +  ageee  ®  ee  +  a^z^z  ®  e^;  (4.36) 

Grr  =  ^{(A  +  2^)Uq  +  X{uo/r  +  Wg  +  Sg)} 

+  k^{(X  +  2fj,)(u[  -  v!g  /2)  +  A(ui/r  +  -  \  [{ug/r)'^  +  w'i^  +  sjf]  -  2w'gUg/r) 

+  \{vi+  6i>2  +  ^Vg)u2  +  \{vi+  2v2)[{wg  +  Mg/r)^  +  Sg^ 

+  2Mo(itg/r  +  uig  +  Sg)]  +  i'is'g{uo/r  +  w'g) 

-  [(A  +  2/t)uo  +  X{uo/r  +  Wg  +  Sg)](3Mo  +  uo/r  +  Wg  +  Sg)},  (4.37) 

<J8S  =  k{{X  +  2^){wg  +  Mg/r)  +  X(ug  +  Sg)} 

+  fc^{(A  +  2^)(Mi/r  -  A[(Mg/r)2  +  Wg^]  -  2w'gUg/r)  +  A(m'i  -  u'i} /2  +  sj  -  SoV2) 

+  |(z>i  +  6t>2  +  8i>3)(M;o  +  Mg/r)^ 

+  |(i>i  +  2t>2)[uo^  +  s'^  +  2{u'g  +  So)(wo  +  ug/r)]  +  hu'gs'g 

-  [(A  +  2^i){wg  +  Ug/r)  +  A(mo  +  So)](wo  +  3ug/r  +  Swg  +  Sg)},  (4.38) 

(Zzz  =  fc{(A  +  2/r)sg  +  A(uo  +  Wg  +  ug/r)} 

+  {{X  +  2fi){s[  -  Sg/2)  +  X{u[  +  ui/r  -  +  {ug/rf  +  w'^]  -  2w'gUg/r) 

+  \{vi+  6i'2  +  8t>3)so2 

+  A(i>i  +  2£'2)[mo^  +  {w'g  +  Ug/rf  +  2{w'g  +  Mg  +  Mg/r)so]  +  V]_u'g{w'g  +  Ug/r) 

-  [(A  +  2/r)so  +  A(mo  +  ug / r  +  w'g)\{u'g  +  ug/r  +  w'g  +  Ssg)}-  (4.39) 
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These  stress  components  are  consistent  with  (3.36).  Applying  (A. 39)  in  spatial  physical  components,  equilibrium  equations 
V  •  (T  =  0  reduce  to 


drCJrr  +  (crr  “  (T99)/r  =  0,  deage  =  0,  d^a^^  =  0.  (4.40) 

Since  tCg  =  constant,  Cauchy  stress  depends  potentially  only  on  r  and  z,  so  the  second  of  (4.40)  is  trivially  satisfied. 

First  consider  equilibrium  conditions  corresponding  to  stress  components  linear  in  k,  i.e.,  the  first  of  (3.37),  where 

aoi'r,z)  =  [(A  +  2^)uo  +  X{uo/r  +  Wg  +  So)]er  0 

+  [(A  +  2^i){w'q  +  utg/r)  +  A(ug  +  So)]^^  0  eg 

+  [(A  +  2/r)sQ  +  A('Ug  +  uq/t  +  '(Ug)]e2  0  e^.  (4.41) 

The  first  order  (i.e.,  linear  in  k)  axial  equilibrium  equation  produces,  with  so(0)  =  0  to  fix  translational  invariance, 

dz{aQ)zz  =  0  ^  Sg  =  0  ^  So  =  C3Z,  (4.42) 

where  cg  is  a  constant  determined  by  boundary  conditions  on  ends  of  the  cylinder.  Thus  (To  =  The  first  order  radial 

equilibrium  equation  becomes 

(A  +  2/r)(uo  +  Uq/t  -  uo/r"^)  =  2fj.WQlr.  (4.43) 


This  is  an  inhomogeneous  second-order  ordinary  differential  equation  of  Cauchy-Euler  type,  with  general  solution  to  the 
homogeneous  equation  of  the  form  ci  r  +  C2 /r  and  total  solution 


uo  =  Cir  +  cg/r  +  [^w'q/{X  +  2^)]r  Inr. 

Applying  boundary  conditions  (crg)rr|r=ro  =  ~P  (o-gj^r  |r=flo  =  0^  the  solution  is 

,,  r2 

Cl  =  x(l  -  2(2)  -  w'q/2  -  VC3  - 


^  ,  Rq  In  Ro  -  rg  In  rg 

^0  R2  _  2 

-^0 


A  “h  2ifi 


'0 


C2  =  xRl  -  X  =  b/(2M)][co/(^o  -  ’’o)]- 

A  -h  Z/i,  /Xq  f  q  IQ 


Substituting  into  (4.44),  radial  displacement  in  the  linear  solution  is 


(4.44) 

(4.45) 

(4.46) 


Uo  = 


(u  w  1  -  2(2  lni?o  -  rg  Inrg 

- 1"  (1  -  2z2)x  -  VC3  -  - - - - - K - 

47r  47r  1  —  z/  R?.  —  r? 


X-^0  1  ,,  d2  „2 


"0  '  0 
uj  1  —  2v 
47r  1  —  12 


'  In  r. 


Constant  cg  =  sg  can  now  be  determined  for  either  of  the  conditions  in  (4.23).  For  vanishing  average  axial  stress, 

p9  =  +  7r  pr  =  Ro 


pa  =  -tn  pr 

J  0=  —71-  J  r  = 


{ao)zzrdrdO  =  0  cg  =  -2ux/il  +  v), 


(4.47) 


(4.48) 


meaning  that  positive  core  pressure  leads  to  axial  contraction  for  Z2  >  0.  For  the  infinitely  long  disclination  line,  cg  =  0, 
and  the  linear  elastic  stress  field  is 


(1  w  fj.  UJ  Rl  In  Ro-rl  In  ro 

(cro)rr  =  j - - ^2 - 2 - 

1  —  z/  ZTT  1  —  Z^  ZTT  fCn  —  Tn 


^0 


w  i?grg  Ro 


H)  -  '  g  ro 

32 


1 

^  +  2(iX) 


fi  UJ  n  w  i?g  In  i?g  -  rg  In  rg 

{cTo)99  =  z. - —  (lnr+  1)  -  - - - - -2 - 2 - 

1  —  (2  27r  1  —  (2  27r  Rf-.  — 


(4.49) 


+ 


M  ^  ^ocg 
2mx«o  -  ) - - 


1 

1  -  (2  27r  i?g  -  rg  rg 


1 

~~2  d~  2^X5 


(4.50) 
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{O'o)zz 


VfJ,  iO 

1  —  v  2tt 


(21nr  +  1) 


21/ fj,  oj  Rq  In  Ro  —  f 0  111  I’D 
1  —  27r  ^0  — 


+  4iy/j,x- 


(4.51) 


This  solution  agrees  with  [30]  when  the  core  is  neglected  such  that  ro  ^  0  and  x  =  0,  in  which  case  strain  energy  per  unit 
length  is  /x(fcaji?o)^/[167r(l  —  v)].  All  components  of  cto  are  linear  in  w  or  x  such  that  k  factors  out  of  the  final  solution 
kao. 

Henceforth,  only  the  infinitely  extended  disclination  line  geometry  is  addressed,  which  corresponds  to  sq  =  si  =0  and 
z  =  Z,  and  the  solution  is  independent  of  z.  First-order  stresses  (linear  in  k)  for  this  case  are  (4.49)-(4.51).  Now  consider 
equilibrium  equations  for  stress  components  quadratic  in  k,  i.e.,  the  second  of  (3.37).  From  (4.37)-(4.39), 


+  2/i)u'i  +  Xui/r]er  0  +  [(A  +  2^)ui/r  +  \ui]ee  0  eg 

-f  [A(ui  +  ui/r)]ez  0  e^, 


(4.52) 


T(r)  =  {  I  (f'l  -I-  6i>2  +  8z>3  -  7A  -  14/x)uo^  +  (Di  +  2i>2  -  4A  -  2fi)uQUo/r 
+  \  {vi+  2z>2  -  3A)(Mo/r)^  +  \  {v\  +  2v2  -  3A)wq^ 

+  (i>i  +  2j>2  —  4A  —  2^)uqWq  +  (i>i  +  2i>2  —  4A)woMo/r}e,.  0 
+  {\{vi+  2C>2  -  3X)uq  +  (z>i  +  2i>2  -  4A  -  2fj,)uQUo/r 
+  \  {vi+  6i>2  +  8i>3  —  IX  —  14/x)(uo/r)^  +  \  {v\  +  6i>2  +  8^’3  —IX—  14^)wq 
+  (Di  +  2v2  -  4A  -  2^)u'qw'q  +  (z>i  +  6i>2  +  8;>3  -  8A  -  16p)woMo/r}ee  0  eg 
“I"  { 21^2  ~  3A)uq^  +  (f'l  —  2X)uqUo /r  +  |  (Di  +  21/2  —  A)  {u^/r 
+  A(i>i  +  2v2  -  X)w'^  +  (j>i  -  2X)uqWq  +  (Di  +  2X)wQUo/r}ez  0  e^.  (4.53) 

The  radial  component  of  second-order  equilibrium  equation  in  (3.37)  is 

dr{(Tl  +  T)rr  +  [(cti  +  T)rr  “  (cti  +  T)ge]/r  =  0.  (4.54) 

Substituting  from  (4.52)  and  (4.53),  this  becomes  an  inhomogeneous  second-order  ordinary  differential  equation  of  Cauchy- 
Euler  type: 


u'l  +  u[/r  -  ui/r^  =  Bi/r  +  B2/r^  +  B^/r^  +  i74(lnr)/r; 

Bi  =  -  — [2ciC4(2i>i  -I-  10i>2  +  8i>3  -  llA  -  16^) 

A  H”  Z^i 

-\-  (^{Zui  -\-  10^'2  “1“  —  13A  —  23/i)  —  {Zi)2  ~\~  —  2A  —  7/i) 

—  4ciu>g(;>2  +  2i>3  —  A  —  4^)  +  2c4w'g(i>i  +  2i>2  —  4A  —  p)], 

4 


Bo  =  - 


B4  =  - 


A  -f  2/i 
4 

X  +  2/J.  ^ 
4 


C2[C4^  -  Wg(i>2  +  21)3  -  X  -  3/i)], 


B3  =  -  '  cl(X  +  5/j.  -  2i>2  -4i>3), 


-C4 


[C4(f’i  -f  4z>2  -f  4i>3  -  ^  A  -  8/x)  -  Wg(i>2  +  2i>3  -  X  -  4/4), 


A  -f  2/i 

/i  w  1  —  2i/  oj 


Wg  = 


OJ 


(4.55) 


(4.56) 

(4.57) 

(4.58) 

(4.59) 

(4.60) 


^  /i  -f  2A  27r  1  —  1/  47r  ’  27r 

The  homogeneous  general  solution  to  this  equation  is  of  the  form  Cir  +  C'2/r.  The  particular  solution  is  found  using 


Lagrange’s  method,  which  when  added  to  the  homogeneous  solution  gives  the  total  general  solution 


C  A 

Ui{r)  =  CiT  -\ - -  -\ - ^4-  A2r-lnr  -I-  Agrlv?'  r  +  A4 


In  r 


Ai=B3/8,  A2  =  (2Bi -B4)/4,  A3=B4/4,  A4  = -B2/2. 


(4.61) 

(4.62) 
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Consistent  with  boundary  conditions  and  the  second  of  (4.24),  applying  the  vanishing  second-order  traction  end  conditions 
((Ti)rr  =  —Trr  at  r  =  tq  and  r  =  Rq  gives 


-  ^  +  2(iT7o(r-R!7«‘“" 


+  Di  [(Rl/rl)  In  i?o  -  In  tq]  +  £>5  [(i?o/^o)  Ro  -  In^  Tq] 

-h  De[(lnRo)/rl  -  (lnro)/r^]},  (4.63) 

+  Dsllii*  ft,  -  In’  r„]  +  ft,[(lnft,)/ftj  -  (lnr„)/rj]};  (4.64) 

Di  =  (A  -f  2fijA.2  4“  {i^i  4“  21)2  —  4A)(2ci  4-  —  2/r(ci  4- 

4-  |(i>i  4-  2j>2  —  3X)wq  +  (2i>i  4-  6j>2  4-  4i>3  —  9A  —  9n)cf 

+  (2f'i  4-  8i>2  4-  8O3  —  llA  —  16/r)ciC4  4-  |(£'i  4-  6j>2  4-  8O3  —  7A  —  14p)c4,  (4.65) 

£^2  =  (A  4-  2^)444  4-  2/iC2Wo  —  2(2z>2  4-  4i>3  —  2A  —  7/x)ciC2 

-  (4i>2  4- 8O3  -  3A  -  12^)c2C4,  (4.66) 

D3  =  -2(A  -h  3p)44i  -4  i  (4t>2  +  8z>3  -  2A  -  10^)4 ,  (4.67) 

D4  =  (A  4“  2p)(442  4“  2^43)  4“  A442  4“  2(2)24  4-  6)22  4-  4)23  —  9A  —  9p)ciC4 

4-  2(!2i  4-  2i>2  —  4A  —  4)C4U)o  +  I  (S^'i  4-  181)2  +  16)23  —  25A  —  32fi)cl,  (4.68) 

£>5  =  2(A  4-  p)443  4-  \  {3v\  +  10)22  4-  8)23  —  15A  —  18p)c4,  (4.69) 

Dq  =  — 2/in44  —  2(2)22  4-  4)23  —  2A  —  7/j,)c2C4.  (4.70) 

Second-order  contribution  to  radial  displacement  ui  is  now  fully  determined.  The  total  displacement  and  stress  helds  are 

u{r)  =  kuo{r)  +  k‘^ui{r),  a{r)  =  kao{r)  +  k'^[ai{r)  +  T{r)],  (4.71) 


where  uq  is  given  by  (4.47)  with  C3  =  0,  ui  by  (4.61),  ag  by  (4.49)-(4.51),  cri  by  (4.52),  and  t  by  (4.53).  Second-order 
contributions  are  proportional  to  x^,  or  0;%;  recalling  that  x  c)c  p  =  p/k  and  uj  =  oj/k,  k  does  not  affect  the  hnal 
solution. 

The  complete  solution  depends  on  9  parameters:  geometric  constants  (lj,  rg,Rg),  elastic  constants  (A,  p,  C'i,C'2,  vg),  and 
core  boundary  condition  p.  General  closed  form  expressions  for  stresses  and  displacements  in  terms  of  these  parameters 
have  not  been  obtained  due  to  the  algebraic  complexity  of  the  solution,  but  numerical  values  could  be  found  using  a 
computer  program  for  a  hxed  set  of  parameters.  Energy  per  unit  length  4)  =  27r  similarly  has  not 

been  obtained  in  closed  form,  but  could  be  integrated  numerically.  Further  examination  of  the  solution  using  numerical 
methods  is  reserved  for  future  work. 


5  Point  defect 

5.1  Geometric  description 

Consider  a  deformed  body  05  with  boundary  chB.  In  the  reference  conhguration,  the  image  of  this  body  and  its  boundary  are 
05o  and  c£Bo.  Spherical  coordinates  covering  05  and  05o  are  (r,  6,  p)  and  (i?,  0,  T*)  and  are  described  fully  in  Appendix  B, 
including  basis  vectors;  Cartesian  coordinates  are  (x,  y,  z)  and  {X,  Y,  Z). 

An  Eulerian  description  of  deformation  is  invoked,  in  contrast  to  previous  works  that  used  Lagrangian  descriptions 
[57, 62, 72],  and  a  different  strain  energy  potential  (W  rather  than  W)  is  also  used  here.  The  point  defect  is  treated  as  a 
singular  center  of  dilatation  or  contraction  at  point  p  G  05  located  at  the  origin  of  the  spatial  coordinate  system.  This  defect 
could  physically  correspond  to  a  vacancy,  interstitial,  or  substitutional  atom,  or  an  inclusion  of  small  size.  The  general 
ansatz  for  deformation  due  to  a  point  defect  in  an  isotropic  body  is  of  the  spherically  symmetric  form 

R{r)=r  —  u{r),  Q{0)  =  9,  <i)((^)  =  0.  (5.1) 
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Radial  displacement  is  u(r),  positive  for  dilatation  and  negative  for  contraction.  Rigid  body  translation  is  excluded,  restrict¬ 


ing  m(0)  =  0  such  that  po>  the  referential  image  of  p,  also  occupies  the  origin  of  the  coordinate  system  on  IBq. 

In  curvilinear  coordinates,  the  inverse  deformation  gradient  from  (5.1)  is 

=  F-\^Gr  0  9’-+  F-i®Ge  0  g'  +  0 

=  [1  -  u'{r)]GR  0  s’-  +  Ge  0  g"'  +  0  (5.2) 

Volume  ratio  is 

J-i  =  (Ryr^)  det  F-^  =  (1  -  m')(1  -  u/rf  >  0,  (5.3) 

noting  that  y^Gjg  =  /r^  =  (1  —  ujr)’^ .  Applying  the  first  of  (2.49)  in  curvilinear  coordinates, 

F'l  =  =  F^-\^F^-}^Ga  0  g“-  (5.4) 

The  second  of  (2.50)  applied  to  the  singular  point  defect  is,  with  the  shifter, 

^  ^  F^-\^  =  (5.5) 

=  1  —  5v5{x)5{y)5{z)  =  1  —  5v5{r).  (5.6) 

Inelastic  volume  change  induced  by  the  defect  is  measured  by  scalar  constant  5v  that  can  be  positive  or  negative.  The 
spherical  Dirac  delta  function  is  5{r)  =  5{x)5{y)5{z),  differing  from  that  used  for  cylindrical  coordinates  in  Sect.  3  and 
Sect.  4.  Therefore, 

F{X)^F^{X)  VXgQ3o\Po,  F~\x)  ^  F^-\x)  V®  g  ®\p.  (5.7) 


Vector 'd  entering  Burgers  vector  integral  (2.52)  is  of  the  spherically  symmetric  form  =  L{r)5{r)gr .  For  an  area  element 
of  a  spherical  shell  with  unit  normal  Qr,  the  differential  form  of  (2.52)  is  db(r)  =  i,{r)6{r)rgrda  =  i,{r)6{r)r^  sin  OdOdcfigr. 
Integral  f  db  then  vanishes  identically  for  a  point  defect  at  r  =  0,  but  does  not  always  vanish  for  a  different  coordinate 
system,  similarly  to  the  case  for  the  disclination  in  Sect.  4. 1 . 

5.2  Nonlinear  elastic  analysis  and  general  solution 

A  nonlinear  elastic  boundary  value  problem  for  the  point  defect  is  constructed  as  follows.  Let  body  05  consist  of  an  elastic 
sphere  of  outer  radius  Rq  with  the  ball  0  <  r  <  tq  removed.  Body  5B  is  simply  connected,  as  is  its  image  in  the  reference 
configuration,  So  -  Region  r  <  tq  is  called  the  core  of  the  point  defect,  and  is  elastically  rigid,  i.e.,  F^  =  1  for  r  <  rg. 
Deformation  of  the  core  region  is  inelastic  or  residual: 

j^  =  detF\  (5.8) 

where  the  determinant  applies  in  a  coincident  spatial  basis  {fla  j  g“}  [15].  Spatial  volume  of  the  defect  core  is  no  =  |7rrg ; 
referential  volume  of  the  core  is,  from  (5.6), 

Vo  =  f  J^~^dv  =  f  [1  —  (5n(5(r)]dn  =  no  —  Sv.  (5.9) 

Jvo  Jvo 

Displacement  boundary  conditions  are  imposed  on  the  elastic  body  at  its  inner  surface  r  =  rg  such  that  the  volume  of  a 
spherical  shell  of  radius  rg  and  constant  thickness  i)  is  equal  to  the  volume  change  induced  by  the  core: 

6v  =  dTrrgfj,  f]  =  u{rg).  (5.10) 

Condition  r)  >  0  corresponds  to  a  referential  core  of  radius  smaller  than  ro,  leading  to  radial  expansion  of  the  surrounding 
elastic  medium  in  the  current  configuration.  Conversely,  fj  <  0  corresponds  to  a  referential  core  of  radius  larger  than  rg , 
leading  to  radial  contraction  of  the  surrounding  elastic  medium  in  the  current  configuration.  Requiring  Vq  >  0  leads  to 
constraint  fj  <  rg/3. 

The  solution  procedure  parallels  that  of  Sect.  4.2,  with  spherical  rather  than  cylindrical  symmetry.  Equilibrium  equations 
are  solved  on  05,  which  excludes  region  r  <  rg  and  therefore  excludes  singular  point  p.  From  (5.7),  in  this  domain  the  strain 
energy  function  in  (2.53)  becomes 

W  =  W{F)  =  W[D{F)]=W[E{F)],  (5.11) 
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where  W  is  used  in  the  solution  that  follows.  Equations  (3.28)-(3.31)  apply.  Internal  and  external  boundary  conditions  are 

tt|r  — ro  —  i);  ^\r  —  Rfj  —  \r  —  RQ^Qr  —  0-  (5.12) 

The  first  of  (5.12)  prescribes  radial  displacement  at  the  core  consistent  with  inelastic  volume  change  of  the  core.  The  second 
of  (5.12)  enforces  vanishing  traction  on  the  external  part  of  OIB. 

An  iterative  solution  procedure  in  Eulerian  coordinates  is  invoked.  Equations  (3.34),  (3.36),  and  (3.37)  apply.  Internal 
boundary  conditions  are  of  the  form 

fj  =  kr],  (5.13) 


Displacement  (3.34)  applied  to  (5.1)  results  in 

u{r)  =  kuo{r)  +  k‘^ui{r).  (5.14) 

In  physical  components,  i.e.,  orthonormal  bases  {sa,  Ea},  (5.2)  and  (B.36)  become 


Substituting  from  (5.14), 

{F-\Aa)]  = 


R' 

0 

0 

1-m' 

0 

0 

0 

R/r 

0 

= 

0 

1  —  u/r 

0 

(5.15) 

_0 

0 

R/r_ 

0 

0 

1  —  u/r 

'l  - 

kuQ  — 

0 

0 

0 

1  - 

kuo/r - 

-  k'^ui/r 

0 

,  (5.16) 

0 

0 

1  —  kuQ /r  —  k'^ui/r 

and  (5.3)  becomes 

J~^  =  =  (1  —  kuQ  —  k^u{){l  —  kuo/r—  «  1  —  fc(itQ  +2ito/r).  (5.17) 

In  physical  spherical  components,  strain  tensor  D  is,  from  (3.28)  and  (5.16), 


D  —  Djii^Eji  0  E}i  +  DqqEq  0  Eq  +  0  E^-,  (5.18) 

Drr  =  ku'„  +  k'^{u[-  Dqq  =  =  kuo/r  +  k'^[ui/r  -  ^{uo/rf],  (5.19) 

where  terms  0{k^)  have  been  truncated.  Strain  energy  potential  per  unit  reference  volume  (C.7)  is  used,  truncated  at  third 

order,  leading  to  (3.47)  which  also  applies  here.  Stress  S  =  dW /dD  is,  omitting  terms  of  0{k^), 

S  =  SjfjiEf{  0  Ej^  +  SqqEq  0  Eq  +  S^^E^,  0  E^]  (5.20) 

Srr  =  k{{\  +  2/i)ito  +  2Aito/r} 

+  k‘^{{X  +  2fi){u[  -  /2)  +  2\[ui/r  -  liua/rf] 

+  \{vi+  Qv2  +  +  2(j>i  +  V2){uQlrf  +  2(£'i  +  2v2)u'qUq  /  r} ,  (5.21) 

Sqq  =  =  fc{(2A  +  2^)moA  + 

+  k^{{2\  +  2^)[ui/r  -  \  {uQlrf]  +  \{u[  -  Uq /2) 

+  \{vi+  ‘2v2)uq  +  2(z>i  +  3j>2  +  2i>3)(MoA)^  +  2(z>i  +  z>2)uoMoA}-  (5.22) 

Using  (5.16)  and  (5.17),  Cauchy  stress  a  =  SF~^  is,  omitting  terms  of  0{k^), 

(7  —  OqqGq  Gq  “h  O^^G^  ^(j)t  (5.23) 

(7^^  =  k{{X  +  2fi)uQ  +  2XuQ/r} 

+  A{(A  +  2A(A  -  Mo^/2)  +  2A[ui/r  -  ^{uo/rf] 

+  +  8A) +  2(£'i  +  £'2)(^^o +  2(i)i  +  2£'2)uqUq / r 

—  [(A  +  2/x)uo  +  2Xuo /r]{3uQ  +  2moA)},  (5.24) 
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cree  =  cr^4,  =  k{(2X  +  2fj.)uo/r  +  Xuq} 

+  k‘^{(2X  +  2fi)[ui/r  -  i(Mo/r)^]  +  A(mi  -  u'^ /2) 

+  |(z>i  +  2v2)uq  +  2{vi  +  3z>2  +  2v^){u()/rY  +  2(t>i  +  i'2)uQUo/r 

-  [{2X  +  2fi)uo/r  +  Xuq]{uq  +4moA)}-  (5.25) 

These  stress  components  are  consistent  with  (3.36).  Applying  (B.38)  in  spatial  physical  components,  equilibrium  equations 
W  ■  a  —  0  reduce  to,  in  the  absence  of  shear  stress. 


drarr  +  {2arr  -  <jgg  -  a^^)/r  =  0,  dgagg  +  {agg  -  cot  0  =  0,  d^a^^  =  0.  (5.26) 

Noting  that  a  =  a(r)  and  agg  =  the  second  and  third  equilibrium  equations  are  trivially  satisfied,  and  radial  equilib¬ 
rium  becomes 

darr /(ir  +  2{arr  —  (Jgg)/r  =  0.  (5.27) 

First  consider  equilibrium  conditions  corresponding  to  stress  components  linear  in  k,  i.e.,  the  first  of  (3.37),  where 

o’o(t’)  =  [(A  +  2fj,)uQ  +  2Xuo/r]er  0  +  [(2A  +  2fj,)uo/r  +  Auo](ee  0  eg  +  (g)  e^).  (5.28) 

First  order  radial  equilibrium  equation  (5.27)  becomes 

(A  +  2/r)(u"  +  2u'Jr  -  2uo/r^)  =  0.  (5.29) 

This  is  a  homogeneous  second-order  ordinary  differential  equation  of  Cauchy-Euler  type  with  general  solution 

uo  =  cir  +  C2/r'^ .  (5.30) 


Applying  boundary  conditions  mo(?'o)  =  V  (cro)rr |r=i?o  =  ^  consistently  with  (5.12)  and  (5.13),  the  solution  is 

r  (3A-f2p)rg?? 

^  [(3A  +  2p)  +  4/r(ro/i?o)^]-Ro  ’  ^  3A  +  2/r -f  4p(ro/i?o)^  ’ 

First  order  stresses  become 

(CTo)rr  =  (3A  -f  2p)ci  -  4fj,C2 / ,  {cro)gg  =  =  (3A  -f  2p)ci  -f  2/rc2/r^.  (5.32) 

First-order  pressure  is  —  |tr(To  =  —3Kci,  with  AT  =  A  -f  the  bulk  modulus.  Displacement  uq  and  stresses  (5.32)  are 
linear  in  ry  =  fj/k,  so  k  does  not  affect  the  solution,  which  agrees  with  previous  derivations  [14, 81].  Strain  energy  density 
is 

W=  ifcVo  :Do+0(k^)=  lk"^l(3X  +  2fx)cf  +  4fxcl/r<^]  +  O(k^).  (5.33) 

Integrating  over  iB  and  dividing  by  total  volume  of  the  elastic  sphere  minus  core  gives  the  average  strain  energy  density  in 
the  linear  approximation 

— J-iVFdt;=  |fc2[(3A  +  2M)c?+4/ic2/(Aoro)3]  +  0(fc3).  (5.34) 

47r^KQ  —  TqJ 

Now  consider  equilibrium  equations  for  stress  components  quadratic  in  k,  i.e.,  the  second  of  (3.37).  From  (5.24)  and 
(5.25), 

=  [(^  +  2/J,)u'i  +  2Aui/r]e,.  0  Br  +  [(2A  -f  2fj,)ui/r  +  Xu\]{eg  0eg  +  e^0  e^),  (5.35) 

T(r)  =  {\{vi+  6^2  -I-  8i>3  -IX-  +  2{vi  +  2v2  -  4A  -  2fj.)uQUo/r 

+  (2Fi  -I-  2i>2  -  5X){uo/r)^}er  0 
+  +  2v2  —  3A)uo^  -I-  2(j>i  -I-  j>2  —  3A  —  fi)uQUo/r 

+  (2z>i  -I-  6f’2  -I-  Avs  -  9A  -  9fi){uo/ry}{eg  0  eg  +  0  e^).  (5.36) 
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Second-order  equilibrium  equation  (5.27)  is 

d(cri  -b  T)„/dr  -b  2[(cri  -b  T)rr  “  (ci  +  T)e0]/r  =  0.  (5.37) 

Substituting  from  (5.35)  and  (5.36),  this  becomes  an  inhomogeneous  second-order  ordinary  differential  equation  of  Cauchy- 
Euler  type: 


18(202  +  4O3  —  A  —  5/r)c| 
(A  -b  2^)r^ 


(5.38) 


The  general  solution  to  the  homogeneous  equation  is  of  the  form  Cir  +  C'2/r^;  adding  this  to  the  particular  solution  gives 
the  total  solution 


=C'ir  +  C'2/r2+bli/r5,  Ai  =  [(2z>2  +  dOg  -  A  -  5^)c^]/(A  +  2^).  (5.39) 

Consistent  with  boundary  conditions  (5.12)  and  (5.13),  applying  second-order  end  conditions  tti(ro)  =  0  and  {ai)rr  = 
—Trr  at  r  =  i?o  gives 


Or  — - g- 

’’5 


C2  =  - 


1 


3i^  +  4/r(ro/i?o)3 
1 


iKAi 


Ih 


3K  +  Mro/Rof 


3KAi  „  ,  B-zrf-, 


+  Bir^o  + 


mo 

R3 


El 
Rt 

Rt  J 


Bi  =  —  lcf(90i  +  I8O2  +  8O3  —  33  A  —  22/i), 

B2  =  4ciC2  (3O2  +  4O3  -  3A  -  8^), 

Bo  =  — Cg  [6O2  +  I6O3  —  3A(1  +  Ai  /c2 )  —  10/r(2  +  Ai  /c^ )] . 


(5.40) 

(5.41) 

(5.42) 

(5.43) 

(5.44) 


Second-order  contribution  to  radial  displacement  ui  is  now  fully  determined.  The  total  displacement  and  stress  helds  are 


u{r)  =  kuo{r)  +  k^ui{r),  (j{r)  =  kao{r)  +  k'^[ai{r)  -bT(r)], 


(5.45) 


where  uq  is  given  by  (5.30),  ui  by  (5.39),  cto  by  (5.32),  Ui  by  (5.35),  and  r  by  (5.36).  Second-order  contributions  are 
proportional  to  recalling  that  rj  =  fj/k,  k  does  not  affect  the  final  solution. 


5.3  Solution  for  ideal  crystal  with  Cauchy  symmetry 

The  complete  solution  depends  on  8  parameters:  geometric  constants  (rg,  i?o),  elastic  constants  (A,  /r,  Di,  fzg,  f'3),  and  core 
boundary  condition  Sv  =  47rrgf}.  The  solution  is  further  examined  for  a  reduced  set  of  parameters.  The  geometry  and 
displaced  core  volume  can  he  expressed  in  terms  of  dimensionless  parameters  A  and  rjo : 

A  =  Ro/ro  >1,  rjo  =  rj/rg  =  5v/{i'Krl)  <  1/3.  (5.46) 

Imposing  Cauchy  symmetry  on  the  elastic  constants  as  in  Sect.  3.3,  (3.85)  applies,  which  reduces  independent  elastic 
constants  to  K  and  K' .  The  reduced  set  of  5  parameters  entering  the  solution  is  now  (A,  rg,  ryg,  AT,  AT').  Noting  that  all 
coefficients  of  displacement  depend  only  on  dimensionless  ratios  of  elastic  constants,  normalized  radial  displacement  u/tq 
versus  r/rg  depends  only  the  set  of  3  parameters  (A,  r]o ,  K').  Since  stress  components  are  linearly  proportional  to  an  elastic 
constant,  dimensionless  stress  a  jK  versus  r/rg  also  depends  only  on  this  set,  as  does  normalized  volume  change 

^  ^  (4/3)^A^  -  (4/3)^[Ag  -  u{Ro)f  ^  A^  -  [A  -  rr(Ag)/rg]3 
8v  dTrrgf)  3r]o 

Normalized  radial  displacement  rt/rg  is  shown  versus  Eulerian  radial  coordinate  r  in  Eig.  4.  The  effect  of  core  displace¬ 
ment  is  shown  in  Fig.  4(a)  for  an  elastic  sphere  of  size  A  =  i?g  /rg  =  100  and  K'  =  4:  dilatation  (i.e.,  radial  expansion) 
increases  with  increasing  ryg  but  decreases  rapidly  with  increasing  r.  Recall  that  rjo  >  0  corresponds  physically  to  a  large 
substitutional  atom  or  interstitial,  while  r]o  <  0  corresponds  to  a  small  substitutional  atom  or  vacancy.  The  effect  of  non¬ 
linear  elastic  constant  K'  is  shown  in  Fig.  4(b)  for  rjo  =  0.3  and  rjo  =  —0.3:  dilatation  increases  with  increasing  K'  for 
rjo  >  0  and  radial  contraction  decreases  with  increasing  AT'  for  rjo  <  0.  The  effect  of  sphere  size  i?g  is  shown  in  Fig.  4(c) 
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(c)  (d) 

Fig.  4  Total  radial  displacement  for  point  defect:  (a)  variable  core  displacement  rjo  (b)  variable  K' ,  r/o  =  0.3,  —0.3  (c) 
variable  sphere  size  Rq  =  Aro  with  rjo  =  0.3  (d)  variable  sphere  size  Rq  =  Aro  with  rio  =  —0.3. 


(a)  (b) 

Fig.  5  Cauchy  stress  for  point  defect  with  K'  =  A  and  A  =  100:  (a)  radial  stress  (b)  circumferential  stress. 


and  Fig.  4(d)  for  rjQ  =  0.3  and  r]o  =  —0.3:  for  A  >  100,  radial  displacements  are  indistinguishable.  For  A  <  10,  increased 
magnitudes  of  u  are  observed  since  the  unconstrained  surface  at  Rq  is  closer  to  the  core. 

Normalized  Cauchy  stress  held  a/K  is  shown  in  Fig.  5  for  A  =  100  and  K'  =  4.  Radial  stress  arr  is  shown  in 
Fig.  5(a),  circumferential  agg  —  in  Fig.  5(b).  Maximum  magnitudes  at  r  =  tq  are  on  the  order  of  A'/2.  Radial  and 
circumferential  stress  components  are  of  opposite  sign  and  decay  rapidly  within  r  <  lOro. 

Nonlinear  and  linear  solutions  for  stress  and  normalized  volume  change  are  compared  in  Fig.  6  and  Table  2,  respec¬ 
tively.  Stress  components  shown  in  Fig.  6  correspond  to  hxed  parameters  A  =  100,  rjQ  =  —0.3  (linear  and  nonlinear)  and 
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r/r„ 


K'=8 

K'=4 

- ®rr 

K'=0 

- 

linear 

- *  ®0(l 

,  K’=8 

■  ®(I0 

,  K'=4 

- Oyy 

,  K'=0 

^00 

linear 

Fig.  6  Nonlinear  and  linear  elastic  solutions  for  stress 
field  of  point  defect  with  A  =  100  and  rio  =  —0.3. 


K'  =  0,4,8  (nonlinear).  For  K'  =  4,  stresses  obtained  from  the  nonlinear  solution  are  almost  indistinguishable  from 
corresponding  stresses  of  the  linear  solution.  For  K'  ^  4,  differences  among  stress  predictions  are  evident  near  the  core, 
e.g.,  for  r  <  5ro.  Solutions  for  normalized  volume  change  in  Table  2  are  obtained  using  (5.47).  The  linear  solution  [14, 81] 
is  AV/Sv  =  3(1  —  tz)/(l  +  tz);  this  yields  a  value  of  |  when  z/  =  j.  As  shown  in  Table  2,  volume  change  is  strongly 
affected  hy  elastic  nonlinearity.  For  a  material  with  K'  differing  substantially  from  4,  nonlinear  theory  should  be  used  to 
predict  volume  change  from  point  defects,  even  though  stress  components  may  not  differ  significantly  between  linear  and 
nonlinear  solutions  for  r  >5ro. 


Table  2  Volume  change  AV/Sv  for  point  defect,  A  =  100. 


VO 

K'  =  0 

K'  =  2 

K'  =  4 

K'  =  6 

K'  =  8 

linear 

0.3 

0.00 

0.93 

1.85 

2.78 

3.70 

1.80 

0.1 

1.20 

1.51 

1.82 

2.13 

2.43 

1.80 

-0.1 

2.40 

2.09 

1.78 

1.47 

1.17 

1.80 

-0.3 

3.60 

2.67 

1.75 

0.82 

-0.10 

1.80 

6  Conclusions 

A  nonlinear  differential  geometric  framework  of  defective  crystals  has  been  presented,  incorporating  a  multiplicative  de¬ 
composition  of  the  deformation  gradient  into  up  to  three  terms  and  an  additive  decompostion  of  a  linear  connection  de¬ 
scribing  covariant  derivatives  of  a  field  of  lattice  director  vectors.  This  theoretical  framework  has  been  applied  to  describe 
three  different  singular  defects:  the  screw  dislocation,  the  wedge  disclination,  and  the  point  defect.  Inelastic  deformation 
is  quantified  for  each  defect  using  generalized  functions  in  the  sense  of  Gel’fand  and  Shilov,  wherein  singular  parts  of 
the  deformation  gradient  are  defined  as  inelastic.  This  particular  application  of  the  multiplicative  decomposition  of  the 
deformation  gradient  for  singular  defects  in  the  nonlinear  theory  is  original.  Analytical  second-order  solutions  have  been 
obtained  in  the  context  of  isotropic  compressible  nonlinear  elasticity.  These  solutions  are  all  new,  incorporating  an  elastic 
potential  depending  on  an  Eulerian  strain  measure  in  material  coordinates  with  elastic  constants  up  to  third  order.  The  iter¬ 
ative  solution  procedure  used  herein  omits,  in  the  stress  components,  products  of  orders  three  and  higher  in  certain  strain 
gradient  components.  Thus  the  present  solutions  are  expected  to  be  reasonably  accurate  when  such  higher-order  products 
are  small,  though  exact  solutions  are  not  available  for  comparison.  For  the  screw  dislocation,  radial  displacements  and 
dilatation  are  strongly  affected  by  elastic  nonlinearity  and  core  pressure,  while  stress  fields  tend  to  converge  to  linear  elastic 
solutions  within  a  radial  distance  of  10-20  lattice  parameters  from  the  core.  For  the  wedge  disclination,  effects  of  core 
pressure  on  radial  displacement  and  stress  fields  that  have  not  been  presented  in  previous  linear  or  nonlinear  analyses  are 
included.  For  the  point  defect,  radial  displacement  and  stress  decay  rapidly  with  distance  from  the  core,  but  total  volume 
change  in  the  body  due  to  the  point  defect  depends  strongly  on  elastic  nonlinearity. 
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Appendix  A:  cylindrical  coordinates 

Reference  cylindrical  coordinates  {X^}  are 

{X\X\X^)  =  iR,e,Z).  (A.l) 

Let  X  denote  the  position  vector  measured  from  the  fixed  origin: 

X{R,Q,Z)  =  RGr  + ZGz.  (A.2) 

Natural  basis  vectors  are 


GA=dAX;  GR{e)  =  dRX,  Ge{R,e)  =  deX,  Gz  =  dzX.  (A.3) 

The  metric  tensor  with  components  Gab  =  Ga  ■  Gb  and  its  inverse  are 


\ 

0 

o’ 

\ 

0 

o’ 

[Gab]  = 

0 

0 

,  [G^’']  = 

0 

0 

0 

0 

1 

0 

0 

1 

Contravariant  basis  vectors  are 

G^  =  G^^Gb-,  G^{e)  =  GR,  G®(i?,0)  =  (l/i?2)Ge,  G^  =  Gz.  (A.5) 

Physical  (dimensionless  unit)  basis  vectors  are 

Ea=GaIV^;  Er{Q)  =  Gr,  £;e(0)  =  (l/i?)Ge,  Ez  =  Gz.  (A.6) 

Since  cylindrical  coordinates  are  orthogonal  (Ea  =  E‘'^),  there  is  no  need  to  distinguish  between  contravariant  and  covari¬ 
ant  physical  components.  Cylindrical  coordinates  are  related  to  Cartesian  coordinates  as  follows: 

X  =  Rcose,  y  =  i?sin0;  X  =  XEx +YEy  +  ZEz;  (A.7) 

Er  =  cos  QEx  +sin0£ly,  Eq  =  —  sinQEx  cos  QEy,  Ez  =  Ez-  (A. 8) 

Spatial  cylindrical  coordinate  chart  {a:“}  has  coordinates 

=  {r,9,z).  (A.9) 

Let  X  denote  the  position  vector  measured  from  the  fixed  origin: 

x{r,  0,  z)  =  rgr  +  zg^.  (A.  10) 

Natural  basis  vectors  are 

9a=daX]  gr{0)  =  drX,  gg{r,0)  =  dsx,  g^=d^x.  (A.ll) 


The  metric  tensor  with  components  gab  =  ga  ■  gb  and  its  inverse  are 
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(A.  12) 
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Contravariant  basis  vectors  are 

g^{0)=gr,  g\r,9)  =  {l/r^)ge,  g^=g.. 

Physical  basis  vectors  are 

ea^galVg^',  er{9)^gr,  eg{9)  =  {l/r)gg,  e^=g^. 

Cylindrical  coordinates  are  related  to  Cartesian  coordinates  as 

X  =  r  cos  9,  y  =  r  sin  9;  x  =  xe^  +  yey  +  ze^ ; 

=  cos  9ex  +  sin  ■,  =  ~  sin  9e.j:  +  cos  9ey ,  62=6^. 

In  this  paper,  coincident  Cartesian  frames  are  used  such  that  Eji  =  5\ea'. 

Ex  =  e-xi  Ey  =  By,  Ez  =  62 . 

Mixed  variant  components  of  the  shifter  tensor  between  reference  and  spatial  systems  are  [14,  15,  35] 
g\{x,X)  =  {g\GAY  g’^  {x)  =  g\G^  {X)- 
where  (•,  •)  denotes  a  scalar  product.  Using  (A.17),  it  follows  that 

9e{0,Q)  =  {g'-{9),GR{Q))  =  e,.{9)  ■  Ej,{e) 

=  (cos  062:  +  sin0ej,)  •  {cosQEx  +  sin  QEy)  =  cos{9  —  0). 


Applying  similar  calculations,  the  matrix  of  components  of  g'^  (r,  9;  R,  0)  is 


9r 

ffe 

9z 

cos(0  —  0) 
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The  inverse  of  \g\  ]  is 

[gt] 

,  such  that 

II 

'gf 

9i 

gf 

cos(0  —  0) 
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1 

Letting  G  =  det[GAB]  and  g  =  det]^^;,],  determinants  obey 

det[g^]  =  =  R/r,  det[g^]  =  \fgjG  =  rjR. 


In  cylindrical  coordinates,  deformation  x  =  x(X)  is  of  the  form 

r  =  r(R,0,Z),  9  =  9(R,e,Z),  z  =  z(R,0,Z). 

Referred  to  natural  bases  {ga ,  G'^},  the  deformation  gradient  is 
E  =  dAX^G^=  dAx'^ga  ®G^  =  F^^ga  ® 

=  duT  gr  ®  G^  +  Oqt  gr  Z>  G®  +  dyr  gr  ®  G^ 

+  d]i9  gg®G^  F  3^9  gg®G^+  dz9  gg  <Z  G^ 

+  3rz  ^2  <8  +  9ez  ^2  C)  G®  +  dzz  g^®G^. 

Referential  and  deformed  volume  elements 

dU  =  VCdX^  dX^  dX^  =  RdRd0dZ ,  du  =  ^/gdx^dx^dx^  =  rdrd9dz 


(A.  13) 

(A.  14) 

(A.  15) 
(A.  16) 

(A.17) 

(A.  18) 

(A.  19) 

(A.20) 

(A.21) 

(A.22) 

(A.23) 


(A.24) 

(A.25) 
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are  related  through  the  Jacohian  determinant 

J  =  Av/AV  =  Aet[F'^j^]yJ g / G  =  det[9Aa;“]  det[gf  ]  =  Aei[dAX°']{r / R) .  (A.26) 

Using  (A. 5),  (A. 6),  and  (A.  14)  to  convert  natural  bases  to  orthonormal  bases  {sa ,  Ea  },  the  two-point  deformation  gradient 
in  physical  components  is 

F  =  F(^aA)ea  <^Ea 

=  dj^r  Er  0  E]i  +  (l/i?)(9er  0  Eq  +  dzr  ®  Ez 

+  rduO  eg®  Er  +  {r/R)de9  eg  Eq  +  rdzO  eg  ®  Ez 

+  Brz  e^®ER  +  {l/R)dez  e^®  Eq  +  dzz  e^,®  Ez,  (A.27) 

where  physical  scalar  components  are  written  in  angled  brackets.  The  Jacobian  is 

J  =  det[F(,^)].  (A.28) 

Using  the  shifter  (A.21),  the  deformation  gradient  can  be  expressed  completely  with  respect  to  Lagrangian  bases: 

F  =  dAX^ig^Ge)  ®G^  =  ]Gb  0  G^ 

=  [dRrcos{9  —  0)  —  rd]i9sm{6  —  0)]Gjj  0  G^ 

+  [90r  cos(0  —  0)  —  r9e0sin(0  —  0)]Gij  0  G® 

+  [9^rcos(0  —  0)  —  rdz9sm{9  —  0)]Gij  0  G^ 

+  [{l/R)dBrsm{9  -  0)  -f  {r/R)dR9cos{9  -  0)]Ge  0  G^ 

+  [{1/R)dersm{9  -  0)  +  {r/ R)de9 cos{9  -  0)]Ge  0  G® 

+  [{1/ R)dzrs\n{9  —  0)  -f  {r / R)dz9 cos{9  —  0)]Ge  0  G^ 

+  [da  zj  Gz  0  G^  +  [ae  z]  Gz  0  G^  +  [dz  z]  Gz0G^.  (A.29) 

In  physical  components  this  becomes 

E  =  F(^ab)Ea  0  Er 

=  [9flrcos(0  —  0)  —  rdR9sm{9  —  Q)]Er  0  Er 
+  [(l/i?)9er  cos(0  —  0)  —  {r / R)de9 sm{9  —  Q)]Er^  0  Eq 
+  [dzrcos{9  —  0)  —  rdz9sm(9  —  Q)]Er  0  Ez 
+  [Ort siii{9  —  0)  +  rdR9cos{9  —  0)]£l0  0  Er 
+  [(l/i?)90r  sin(0  —  0)  +  {r/R)de9 cos{9  —  0)]£l0  0  Eq 
+  [dzrs\n{9  —  0)  +  rdz9cos{9  —  Q)]Eq  0  Ez 

+  [dRz]Ez  0  Er  +  [{l/R)dez]Ez  0  Eq  +  [dzz]Ez  0  Ez-  (A.30) 

Similarly,  expressing  E  entirely  with  respect  to  natural  spatial  bases, 

F  =  dAX°-ga  0  {ggg'')  =  [dAX°^gg]ga  0  g'’ 

=  [9/jrcos(0  —  0)  +  (l/i?)90r  sin(0  —  &)]gr  0  g^ 

+  [dR9cos{9  —  0)  +  (l/i?)90  0sin(0  —  Q)]gg  0  g^ 

+  [dRZCOs{9  —  0)  +  (l/i?)90zsin(6>  —  0)]gz  0 g’ 

+  [— r9i;rsin(0  —  0)  +  {r / R)dQr  cos{9  —  0)]gr  0  g^ 

+  [—rdR9 siii{9  —  0)  +  {r / R)dQ9 cos{9  —  Q)]gg  0  g^ 

+  [—rdRZsin{9  —  0)  +  (r / R)dez cos{9  —  Q)]gz  0  g^ 

+  [dzr]gr  0g^  +  [dz9]gg  0  g""  +  [dzz]gz  0  gF  (A.31) 
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In  spatial  physical  components, 

F  =  F(^ab)f^a  ®  Sb 

=  [9/jrcos(0  —  0)  +  {l/R)dersin{0  —  0)]er  0 
+  [rdR0cos{6  —  0)  +  {r / R)dQ9 sm(9  —  0)]e6i  0 
+  [d]iZCOs{9  —  0)  +  (l/i?)9e-zsin(0  —  0)]ez  0 
+  [— i9i^rsin(61  —  0)  +  (l/i?)9er cos(0  —  0)]e,.  0  eg 
+  [—rdR9 siii(9  —  0)  +  {r/R)dQ9 cos{9  —  &)]ee  0  eg 
+  [—d]izs\n{9  —  0)  +  {1  / R)dQ z cos,{9  —  Q)\e^  0  eg 

+  [dzr]er  0  +  [rdz9]eg  0  +  [dzz]e2  0  e^.  (A.32) 

Letting  F  — >  F^ ,  (A.29)  becomes  an  example  of  the  plastic  deformation  gradient  referred  to  coincident  curvilinear  coor¬ 
dinate  frames  in  reference  and  intermediate  configurations,  as  used  in  [14,  15,  18,  75].  With  F  F^,  (A.31)  becomes  an 
example  of  coincident  curvilinear  systems  in  current  and  intermediate  configurations,  as  in  [14,  15]. 

In  cylindrical  coordinates,  inverse  deformation  X  =  X(£c)  is  of  the  form 

R  =  R{r,9,z),  0  =  0(r,6»,z),  Z=Z{r,9,z).  (A.33) 

Referred  to  natural  bases  {Ga  ,9°"},  the  inverse  deformation  gradient  is 

F-i  =  daX  0  =  daX^GA  0  =  F-\iGA  0 

=  drR  GriZ  g''  +  dgR  Gr  +  d^R  Gr  0 

+  drQ  Go  0  g"  +  dgQ  Ge  0  +  5,0  Ge  0  g" 

dr  Z  Gz  0  g^  F  dgZ  Gz  ®g^-\-dzZGz®g^-  (A.34) 

The  inverse  of  the  Jacobian  determinant  J~^  =  1/J  is 

j-i  =  dV/dv  =  det[F-\f  ]i/^  =  det[5<jX^]  det[g^]  =  det[5„A^](i?/r).  (A.35) 

In  orthonormal  bases  {Ea,  Ga},  two-point  inverse  deformation  gradient  is 

0  e, 

=  drR  ErZ)  Gr  +  {l/r)dgR  Er  Z)  Gg  +  dzR  Er  0  e, 

-I-  RdrQ  Eq  0  er  +  {R/r)dg<d  Eq  Z)  Gg  +  RdzQ  Eq  0  e, 

-I-  drZ  Ez  <ZGr  +  {l/r)deZ  Ez  Gg  +  d^Z  Ez  ®  g^,  (A.36) 

with  inverse  Jacobian  determinant 

J-^  =  det[F^2a)\-  (A-37) 

A  spatial  displacement  field  in  physical  cylindrical  coordinates  can  be  constructed  using  the  shifter: 

u^x-  X  =  rgr  +  zg^-  RGr  -  ZGz  =  rgr  -  Rg^ga  +  {z  -  Z)Gz 
=  [r  —  Rcos{9  —  0)]er  -f  [i?sin(6l  —  0)]efl  -f  [z  —  0]e,  =  UrGr  +ugeg  +  u^e^.  (A.38) 

The  divergence  of  a  second-order  tensor  field  A{x)  =  Ai^rib)^a  0  ^b  is,  in  physical  cylindrical  coordinates  [65], 

V  •  A  =  [(l/r)5,.(rA„)  +  {l/r)dgAgr  +  d,A,r  -  {l/r)Agg]er 

+  [{l/r)dr{rArg)  +  {l/r)dgAgg  +  5,A,e  -f  {l/r)Agr]Gg 

+  [(l/r)dr{rArz)  +  {l/r)dgAe^  -f  (A.39) 


(§)  2015  WILEY-VCH  Verlag  GmbH  &  Co.  KGaA,  Weinheim 


www.zamm-journal.org 


ZAMM  ■  Z.  Angew.  Math.  Mech.  (2015)  /  www.zamm-joumal.org 


5 


Appendix  B:  spherical  coordinates 

Reference  spherical  coordinates  {X^}  are 

{X\X\X^)  =  {R,Q,^).  (B.l) 

Let  X  denote  the  position  vector  measured  from  a  fixed  origin: 

X{R,e,^)  =  RGR{e,^).  (B.2) 

Natural  basis  vectors  are 

GA=dAX;  GR{e,<i>)  =  dRX,  G0(R,0,^)  =  d0X,  G$(i?,0,$)  =  (B.3) 

The  metric  tensor  with  components  Gab  =  Ga  ■  Gb  and  its  inverse  are 


'l 

0 

0 

'l 

0 

0 

[Gab]  = 

0 

0 

,  [G^^]  = 

0 

l/i?2 

0 

_0 

0 

i?2  sin^  0 

_0 

0 

sin2  0)_ 

Contravariant  basis  vectors  are 

G‘^=G^^Gb;  G^  =  Gr,  G®  =  (l/i?2)Ge,  =  [1/(7?^  sin^  9)]^$.  (B.5) 

Physical  (dimensionless  unit)  basis  vectors  are 

Ea=GaIV^\  Er  =  Gr,  E0  =  il/R)G0,  EA.  =  [l/{Rsme)]G„.  (B.6) 

Since  spherical  coordinates  are  orthogonal  (Ea  —  E^),  there  is  no  need  to  distinguish  between  contravariant  and  covariant 
physical  components.  Spherical  coordinates  are  related  to  Cartesian  coordinates  as  follows: 

X  =  i?sin0cosT>,  F  =  i?sin  0  sin  T>,  Z  =  Rcos0;  X  =  XEx +YEy  +  ZEz;  (B.7) 

Er  =  sin  0  cos  ^Ex  +  sin  0  sin  ^Ey  +  cos  QEy , 

E0  =  cos  0  cos  ^Ex  +  cos  0  sin  ^Ey  —  sin  SEz , 

£)$  =  —  sin  ^Ex  +  cos  ^Ey .  (B.8) 

Spatial  spherical  coordinate  chart  {x“}  has  coordinates 

=  {r,9,(l)).  (B.9) 

Let  X  denote  the  position  vector  measured  from  a  fixed  origin: 

x{r,e,(j))  =  rgr{0,(j)).  (B.IO) 

Natural  basis  vectors  are 

ga=daX;  gr{9,4>)  =  drX,  gg{r,9,(j))  =  dgx,  g^{r,9,(j))  =  d^x.  (B.ll) 

The  metric  tensor  with  components  gab  =  da  Qb  and  its  inverse  are 

'i  0  0  1  ' 

[gab\=  0  r2  0  ,  [9“']=  0  l/r2  0  •  (B.12) 

0  0  sin^  (/)  0  0  l/(r^  sin^  cj)) 

Contravariant  basis  vectors  are 

9“  =  g°‘'‘gb]  =gr,  g^  =  il/r'^)gg,  g^  =  [l/(r^  sin^  9)]g^.  (B.13) 
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Physical  basis  vectors  are 

ea  =  Oal^/g^;  er  =  Qr,  eg  =  {l/r)gg,  6^  =  [l/(r  sin  6»)]g^ . 
Spherical  coordinates  are  related  to  Cartesian  coordinates  as 


cc  =  rsin0cos(^,  y  =  rsin0sm(j),  z  =  rcos0;  x  =  xe^;  +  yey  +  ze^; 

Br  =  sin  0  cos  4>ej;  +  sin  0  sin  (pey  +  cos  , 
eg  =  cos  0  cos  (jiex  +  cos  0  sin  (j)ey  —  sin  06^ , 
e^  =  —  sin  tpBx  +  cos  4>ey . 

Using  (A.  17),  the  matrix  of  shifter  components  g'\  {r,  0,  (p]  R,  0,  <i>)  is 


9r 

9e 

9h 

9^r 

9e 

9% 

4> 

4> 

4> 

\9r 

90 

9l 

g'jf  =  sin  0  sin  0  cos{(j)  —  <!>)+  cos  0  cos  0, 
gQ  —  i? [sin  0  cos  0  cos((/)  —<!))  —  cos  0  sin  0] , 
g^  =  R  sin  0  sin  0  sin{(j)  —  <!)) , 

9r  =  (l/r)[cos0sin0cos(^  —  <!))—  sin  0  cos  0], 

5e  =  (i?/r)  [cos  0  cos  0cos(())  —  <!>)+  sin  0  sin  0], 

5^  =  (R/r)  cos0sin0sin(())  —  <&), 

9r  —  —  [sin0/(rsin0)]  sin(0  —  <&), 
g^  =  —  [(i?cos0)/(rsin0)]  sin(0  —  <&), 
g^  =  [(i?sin0)/(r  sin0)]  cos(0  —  <i>). 

Letting  G  =  det[G'yiB]  =  sin^  0  and  g  =  det[5a^]  =  sin^  0,  determinants  obey 

det[(7[4]  =  y/ G/g  =  R^  sin0/(r^  sin0),  det[(/f  ]  =  \fgjG  =  sin0 l{Rp'  sin0). 

In  spherical  coordinates,  deformation  x  =  x(X)  is  of  the  form 

r  =  r(i?,0,d>),  0  =  0(i?,0,$),  p  =  p(R,e,<^). 

Referred  to  natural  bases  G  },  the  deformation  gradient  is 

F  =  dAX0G^  =  dAX^ga  0  =  F^^ga  0 

=  Qrt  gr  0  G^  +  Oqt  gr  0  G®  +  0  G^ 

+  dR0  gg0G^+  de0  g0  0G^  +  5$0  gg  0  G® 

+  dB,<j)  g,i,  0  G^  +  dep  g,i,  0  G^  +  d^ip  g^  0  G'^ . 

Referential  and  deformed  volume  elements 

dU  =  VCdX^ dX^ dX^  =  R^  sin  0di?d0d<i),  dt;  =  y^dx^ dx"^ dx^  =  sin0drd0d0 
are  related  through  the  Jacobian  determinant 

J  =  dv/dV  =  det[F'pj)P\-^ g / G  =  det[9Aa;“]  det[gf  ]  =  det[dAX°‘][r'^  sin0/(i?^  sin0)]. 


(B.14) 

(B.15) 


(B.16) 


(B.17) 

(B.18) 

(B.19) 

(B.20) 

(B.21) 

(B.22) 

(B.23) 

(B.24) 

(B.25) 

(B.26) 

(B.27) 

(B.28) 


(B.29) 

(B.30) 

(B.31) 
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Using  (B.5),  (B.6),  and  (B.14)  to  convert  natural  basis  vectors  to  spherical  orthonormal  bases  {ea,EA},  the  two-point 
deformation  gradient  in  physical  components  is 

F  =  F(^aA)ea  ®Ea 

=  dgr  Br  0  Er  +  (l/i?)9er  0  Eq  +  [l/(i?sin 0)]9$r  0  E^ 

+  rdnO  eg  0  Er  +  {r/R)de9  eg  (g)  Eq  +  [r/(i?sin 0)]c)$0  eg  0  E^ 

+  rsinOdiicj)  e<i>  <8  Er  +  {r/R)  sinOdecj)  e^  (g)  Eq  +  [rsin0/(i?sin0)]9$(/)  (g)  £)$.  (B.32) 

Inverse  deformation  X  =  X{x)  is  of  the  form,  in  spherical  coordinates, 

R  =  R{r,e,(j)),  e  =  e{r,0,(j)),  $  =  $(r,6l,</)).  (B.33) 

Referred  to  natural  bases  {Ga  ,  9°  },  the  inverse  deformation  gradient  is 

F-i  =  d,X  0  =  daX^GA  ®g^=  F-\^Ga  ® 

=  drR  Gfl  (g)  g’’  +  dgR  G/j  (g>  g®  +  d^R  Gr  ® 

+  drS  Go  (g)  g"  +  dgO  Ge  0  g®  +  Gq  0  g^ 

+  dr^  G$  0  g^  +  dg^  G$  0  g^  +  G$  (g)  g’^.  (B.34) 

The  inverse  of  the  Jacobian  determinant  J~^  =  1/J  is 

J~^  =  dU/du  =  det[F~\f  ]\/G/g  =  det[9QX^][i?^  sin0/(r^  sin0)].  (B.35) 

With  respect  to  orthonormal  bases  {ca,  Ea},  the  inverse  deformation  gradient  is 

=  drR  Eri^  Br  +  {l/r)dgR  Er  I^eg  +  [l/(r  sin0)]50ii  Er  (g)  e^ 

+  RdrQ  Eq  (g)  Br  +  {R/r)dg<d  Eq  (g  +  [R/ {r  sm  9)]d^Q  Eq  g  e^ 

+  RsinQdr^  E^  g  +  (R/r)  sin099<i)  E^  g  eg 

-I- [i?sin0/(rsin0)]i9,^<i)  £)$  g  e,^.  (B.36) 

A  spatial  displacement  held  in  physical  spherical  coordinates  can  be  constructed  using  the  shifter  of  (B.17): 

u  ^  X  -  X  ^  rgr  -  RGr  =  rg^  -  RgRga 

=  [r  -  RgR]er  +  [-Rrg^j^]eg  +  [-Rr  sm9gff]e^  =  UrBr  +  ugeg  +  u^e^.  (B.37) 

The  divergence  of  a  second-order  tensor  held  A{x)  =  A(^rib)ea  g  in  physical  spherical  coordinates  is  [65] 

V  •  A  =  {{1  / r^)dr{r^ Arr)  +  [l/(rsin6*)]99(sin0A6)r) 

+  [l/(rsin6i)]5,^A,^r  -  {^/r){Agg  +  A^^)}er 

+  {{l/r'^)dr{r'^Arg)  +  [l/(rsin6>)]i99(sin6»A09) 

+  [l/(rsin6i)]5^A,^g  +  {l/r){Agr  -  cot6A^^)}eg 
+  {{l/r'^)dr{r'^  Ar^)  +  {l/r)dg{s\neAg^) 

+  [l/(rsin6*)]5,^A,^^  +  {l/r){A^r  +  A^g))e^.  (B.38) 

Appendix  C:  elastic  potentials 

Let  E  =  WqX  denote  the  deformation  gradient.  Eulerian  strain  D  and  Lagrangian  Green  strain  E,  both  referred  to  material 
coordinates,  are  dehned  as  [16,  28,  29,  85] 

D  =  i(l  -  F-^E^'^),  E  =  ^{F'^F  -  1).  (C.l) 
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For  a  homogeneous  material,  strain  energy  per  unit  reference  volume  is  of  the  form 

W  =  W{F)  =  W[D{F)]  =  W[E{F)].  (C.2) 


First  Piola-Kirchhoff  stress  and  static  linear  momentum  balance  in  the  absence  of  body  force  are 

P  =  dW/dF,  Vo-P  =  0  (VoaP''^=0).  (C.3) 

Thermodynamic  conjugates  to  strains  are  symmetric  stresses  in  material  coordinates 

S  =  dWjdD,  S  =  dWjdE.  (C.4) 

From  symmetries  of  tensors  in  (C.4),  the  chain  rule  and  dF~^j^/dF^g  =  —F~^^F~^^, 


P  = 


dW  dD 

Im  ■  1)F 


Cauchy  stress  then  follows  as  [16,  28,  29,  85] 


^  tW  .dE^ 
dE  '  dF 


a  =  J-^PF'^  =  J-'^F-'^SF~^  =  J-^FSF'^. 


(C.5) 


(C.6) 


Using  Greek  indices  to  denote  Voigt  notation  in  a  Cartesian  frame  (e.g.,  a  =  1,  2, ...  6),  strain  energies  are  written  as 
Taylor  polynomials  in  either  strain  measure: 

W  =  \CafiDaDfj  +  +  •  •  •  ,  (C.7) 

W  =  ^CapEaEfj  +  ^CaPjEaE/jE^  +  •  •  •  .  (C.8) 

Setting  moduli  W /dFdF  =  d^W /dFdF  at  the  unstrained  reference  state,  second-order  elastic  constants  are  equiva¬ 
lent  among  Eulerian  and  Lagrangian  representations: 

Cap  =  Cap  =  Cap.  (C.9) 

Setting  moduli  d^W /dFdF dF  =  d^W /dF dFdF  at  the  unstrained  reference  state,  third-order  elastic  constants  are 
related  as  follows  in  Cartesian  tensor  notation  [16]; 

CijkLMN  =  CijkLMN  +  SjkCjlmN  +  SilCjkMN  +  SimCkLJN  +  SinCkL.JM 
+  ^JkCjlmN  +  SjmCjnKL  +  SjlCjkMN  +  SjnCjmKL 
+  ^KmCijlN  +  SknCjjlM  +  SlmCjjkn  +  SlnCjjkM  ■  (C.IO) 


Relations  identical  to  (C.9)  and  (C.IO)  can  also  be  derived  using  the  expansion  E  =  D  +  2D^  +  ■  ■  ■  and  equating  like 
terms  in  W  and  W  to  third  order  in  strain  components  [66,  90].  For  an  isotropic  solid  [14,  81], 

Cijkl  =  ^5ij5kl  +  K^ikSjl  +  Sil^jk),  (C.ll) 

Cijklmn  =  vi[5ij5kl5mn] 

+  V2[^I.j{5kM^LN  +  SknSlm) 

+  Skl{SimSjn  +  SinSjm)  +  5mn{5ik6jl  + 

+  +  5jmSlm)  +  Sjl{SimSkN  +  ^IN&Km) 

+  Sil{SjmSkn  +  SjnSkm)  +  Sjk{SimSln  +  SinSlm)],  (C.12) 

Cijklmn  =  ^i[5ij6klSmn] 

+  v2[6ij{6kmSln  +  Skn^lm) 

+  Skl{SimSjn  +  SinSjm)  +  Smn{SikSjl  +  SjjSjk)] 

+  vz[Sik{SjmSln  +  SjnSlm)  +  Sjl{SimSkn  +  SinSkm) 

+  Sil{SjmSkn  +  SjnSkm)  +  Sjk{SimSln  +  SinSlm)]-  (C.13) 
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These  also  apply  pointwise  in  a  physical  components  in  an  orthogonal  system,  and  contravariant  components  can  be  recov¬ 
ered  for  general  curvilinear  systems  by  setting  6ij  .  Second-order  constants  obey  the  familiar  relations 


Cii  —  A  +  2/i,  Ci2  —  A,  C44  —  fi]  V  —  A/ (2  A  +  2^) ,  K  —  A  +  ^  /i. 

Third-order  constants  obey  the  following  relations  (CajS-y  or  Cap-y)' 

Cm  =  v\  +  6u2  +  8tz3,  C112  =  vy  +  2v2,  C123  = 

Ci44  =  1^2,  (^155  =  +  2:^3,  C456  =  1^3. 

For  isotropic  materials,  (C.IO)  reduces  to 

t>i  =  Pi ,  P2  =  f'z  +  2A,  Vd  =  i^3  +  3p. 


(C.14) 


(C.15) 


(C.16) 


For  hydrostatic  loading  (cr  =  —pi,  F  =  J^/^1),  pressure  derivative  of  tangent  bulk  modulus  B  at  the  reference  state  is 
[16,  42,  63] 


B'o  =  {dB/dp)\p^o  =  +  2P2  +  |i>3)  +  4  =  +  2123  +  |j^3)- 


(C.17) 


Since  typical  crystalline  solids  have  2  <  B'  <  7  [45],  third-order  constants  are  generally  smaller  in  magnitude 
than  Cap^y,  suggesting  that  series  (C.7)  converges  faster,  i.e.,  requires  fewer  higher-order  terms  for  a  given  accuracy,  than 
(C.8).  Under  hydrostatic  loading,  stress-strain  constitutive  relations  degenerate  to  the  following  pressure-volume  equations 
of  state  [16,  84]: 


p  =  -dW/dJ  =  +  |(B'  -  4)(J-2/3  -  1)],  (Eulerian) 

p  =  -dW/dJ  =  -  |B'(J2/3  -  1)].  (Lagrangian) 


(C.18) 


(C.19) 


Birch-Murnaghan  equation  of  state  (C.18)  is  often  used  to  accurately  fit  hydrostatic  and  shock  compression  data  to  high 
pressures  [11,  45],  while  Lagrangian  equation  of  state  (C.  19)  is  usually  deemed  less  accurate  [45].  For  an  isotropic  material 
also  obeying  Cauchy’s  symmetry  relations  [14,  87], 


X  =  ^  K  v  =  1 


4  ’ 


=  1^3  =  -i^KB'yy. 


(C.20) 

(C.21) 


Applying  (C.16)  with  (C.20)  and  (C.21), 

=  P2  =  -f(|B;-6),  i),  =  -f{lB'o-9).  (C.22) 

It  is  instructive  to  compare  Eulerian  £)-based  theory  with  Murnaghan’s  theory  [61]  based  on  Almansi  strain  e: 

W  =  W[e{F)],  e=  (C.23) 

Because  e  is  referred  to  spatial  coordinates  and  is  not  rotationally  invariant,  strain  energy  W  is  restricted  to  isotropic 
materials,  in  contrast  to  W.  Stresses  are  [61,  84] 


J-^F~'^F-^S,  S  =  dW/de. 

Typically  W  is  expressed  as  a  function  of  three  scalar  invariants  [38,  52,  61]: 


W{Ii,  I2,  I3)  —  {h\  +  fl)Ii  —  2/4/2  +  +  6P2  +  8p3)/f  —  2(^2  +  2p3)/i/2  +  4P3/3; 


h  =  tre, 


I2  =  “  tr(e2)],  I3  =  dete. 


(C.24) 


(C.25) 

(C.26) 


Second-order  constants  A  and  /t  are  identical  to  those  already  introduced.  Third-order  constants  are  usually  related  to 
Lagrangian  constants  by  writing  W  in  terms  of  invariants  of  E,  using  algebraic  relationships  between  these  invariants  and 
those  in  (C.26)  [61],  and  comparing  coefficients  of  like  terms  in  W  and  W,  giving  [38] 


+  44423),  i>2  =  2(A  -  422  -  3423), 


i>3  =  J^3  +  3/i. 


(C.27) 
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From  (C.16)  and  (C.27),  in  general  ^  z>i  and  v-i  ^  t>2,  but  z/3  =  1)3.  These  differences  arise  because  of  different 
truncations  used  to  relate  Eulerian  to  Lagranglan  strain  components  or  invariants.  Note  that  when  ¥  (and  hence  is 
symmetric,  T3  =  e  in  a  common  coordinate  frame,  and  D  and  e  share  the  same  invariants  regardless  of  rotation.  Therefore, 
for  an  isotropic  solid  subjected  to  a  given  F,  setting  Fi  =  i>i,  F2  =  i>2,  and  F3  =  will  produce  W  =  W  and  result 
in  the  same  Cauchy  stress  a  from  (C.24)  or  the  second  of  (C.6),  though  generally  S  ^  S.  A  Birch-Murnaghan  equation 
of  state  also  follows,  i.e.,  p  =  —dW  jdJ  reduces  to  (C.18).  However,  because  their  third-order  constants  are  not  required 
to  be  equal,  in  general,  energy  and  stress  predicted  by  Eulerian  theory  based  on  D  and  that  based  on  e  can  differ,  even 
for  isotropic  materials,  and  only  ID-based  theory  can  be  applied  to  anisotropic  materials,  of  course  removing  symmetry 
restrictions  in  (C.l  1)  and  (C.12). 
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